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SECTION  I 


INTRODUCTION 


1 . OVERV I EW 

The  analysis  of  powered  lifting  systems  of  the  jet  flap  type  requires  an 
accounting  of  the  mutual  aerodynamic  interference  between  the  airframe  (hard 
boundary)  and  the  jet  efflux  (soft  boundary).  This  differs  from  the  classical 
Interference  problem  in  that  boundary  conditions  of  tangential  flow  must 
additionally  be  satisfied  over  a surface  whose  unknown  geometry  is  self- 
adaptive  in  terms  of  both  jet  properties  and  hard  surface  loading  conditions, 
where  such  loading  is  partially  defined  by  the  jet  interference. 

The  development  of  a very  generalized  analysis  technique  is  hardly 
feasible.  However,  the  Externally  Blown  Flap  (EBF),  Upper  Surface  Blowing 
(USB),  Internally  Blown  Flap  (IBF),  and  the  Augmentor  Wing  (AW)  constitute  a 
family  of  powered  lifting  systems  which  are  reducible  to  a common  overall 
aerodynamic  concept,  but  which  differ  primarily  in  their  mechanical  implemen- 
tation and  secondarily  in  aerodynamic  detail.  This  family  Is  best  thought  of 
as  analogous  jet  flapped  wings  differing  primarily  in  the  spanwise  distribu- 
tion of  the  f-fflt  . momenta  at  the  trailing  edge. 

This  report  presents  a relatively  simple  technique  based  on  careful 
application  of  iineari_ed  theory.  The  method  applies  to  the  analogous  systems 
mentioned  ab'-ve.  The  comparison  of  theoretical  results  with  experimental  data 
shows  that  the  method  is  sufficiently  accurate  to  be  very  useful  in  practice. 


2.  JET  FLAP  BACKGROUND 

The  basic  concept  of  our  treatment  is  the  same  as  that  of  a number  of 
other  theoretical  approaches.  I Some  of  their  feature-  will  be  described  sub- 
sequently.) The  jet  Is  considered  as  a thin  sheet  of  fluid  emerging  from  the 
trailing  edge  of  the  wing.  The  flows,  both  internal  and  external  to  this  Jet 
sheet,  are  considered  as  inviscid  so  thf  . no  mixing  or  dlffut'on  takes  place. 
Due  to  the  pressure  differences  exerted  by  the  outer  flow  on  the  jet  sheet  It 
Is  bent  back.  So  far  this  s'.ape  is  unknown.  According  to  this  model  the 
effect  of  the  jet  sheer  on  the  external  field  Is  solely  given  by  its  geometry; 
one  could  think  of  it  as  rigid.  The  pressure  distribution,  which  determines 
Its  shape,  it  determined  by  the  conditions  of  dynamic  equlllorlum  with  the 
outer  field.  Because  of  the  pressure  T’fference  whl^h  exists  between  the  two 
sides  of  the  sheet,  It  can  be  represented  by  a vortex  tayer.  All  tratlirg 
vortices,  those  from  the  wing  as  well  -s  those  from  the  vortex  distribution 
In  the  sheet,  lie  on  this  sheet. 

One  should  remember  that  the  forces  exerted  on  this  sheet  from  the 
outside  are  exactly  balanced  by  the  momentum  cHnge  of  the  jet  within.  The 
forces  exerted  by  the  sheet  on  the  outer  flow  do  not  act  directly  on  the  wing. 
A direct  effect  on  the  wing  is  giver*  by  the  reaction  of  the  jet  as  it  exits 
from  the  wing,  usually  at  its  trailing  edge.  It  is  dete.mined  by  the  momentum 


coefficient  of  the  jet  and  its  direction.  The  pressure  forces  (in  other  words 
the  vortex  distribution  over  the  jet  sheet)  exerted  by  the  jet  sheet  on  the 
outer  flow  cause  a change  of  the  pressure  distribution  over  the  wing.  The 
assessment  of  this  effect  is  the  crucial  step  in  the  treatment  of  jet  flaps. 

It  is  obvious  that  the  vortex  distribution  in  the  sheet  at  stations  close  to 
the  wing  is  of  primary  importance.  In  practical  computations,  one  frequently 
prefers  to  determine  the  combined  lift  of  the  wing  and  of  a suitable  part  of 
the  jet  rather  than  the  lift  on  the  wing  alone.  In  this  case  the  redirected 
jet  momentum  as  it  emerges  from  the  "truncated"  sheet  is  taken  into  account. 
This  merely  amounts  to  a different  kind  of  bookkeeping,  which  is  made  possible 
by  the  fact  that  the  outer  forces  on  the  sheet  are  balanced  by  the  momentum 
changes  of  the  jet.  Sometimes  It  may  be  desirable  to  go  back  to  primary 
quantities,  namely,  the  lift  on  the  wing  by  itself  and  the  momentum  of  the  jet 
as  it  exits  from  the  wing. 

If  one  uses  this  model,  then,  for  an  infinite  swept  wing,  the  jet  sheet 
will  have  the  shape  of  a (general)  cylinder  generated  by  straight  lines 
parallel  to  the  line  along  which  the  Jet  emerges.  Of  course,  for  a wing  of 
finite  dimension  this  configuration  will  appear  only  iocally.  But  this  is  the 
region  in  which  the  vortex  distribution  on  the  jet  sheet  is  of  most  importance. 
Therefore,  one  can  hope  that  by  this  idea  the  effect  of  the  sheet  vorticity  on 
the  wing  is  rather  well  approximated. 

The  classical  mathematical  treatment  of  the  two  dimensional  jet  flap 
problem  is  best  typified  by  that  of  Spence'^,  which  was  corroborated  by 
Malvard's  rheoelectric  anaiopyU),  The  physical  model  consists  of  a wing  with  a 
simple-hinged  mechanical  flap  from  which  a thin  jet  of  fluid  Issues  parallel 
to  the  flap  at  the  trailing  edge.  Flows  both  external  to  and  internal  to  the 
jet  are  invlscld,  so  that  no  mixing  or  diffusion  takes  place.  The  jet  bends 
rearward  due  to  a pressure  differential  across  it  which  results  from  its 
Interference  with  the  external  field.  The  loading  over  the  jet  is  replaced  by 
a distribution  of  vorticity.  The  jet  then  becomes  a vortex  layer  over  which 
the  normal  boundary  conditions  «pply  when  the  jet  and  external  field  are  in 
dynamic  equilibrium. 

The  condition  for  dynamic  equilibrium  is  that  pressure  forces  across  the 
jet  are  balanced  by  centrifugal  forces  on  each  element  of  efflux,  which 
relates  local  vorticity  distribution  to  jet  radius  of  curvature  and  thus  to 
z",  where  z Is  the  distance  from  the  jet  to  the  wing  reference  plane  (small 
slopes  assumed).  This  condition,  together  with  the  normal  thin  airfoil  formu- 
lation, leads  to  a pair  of  Integro-di fferentlal  equations  having  as  unknowns 
(a)  the  vertex  distribution  of  the  airfoil  and  (b)  the  jet  parameter  z"  (and 
hence  the  jet  vortex  distribution).  The  calculated  jet  boundary  Is  required 
to  conform  to  the  injection  angle  at  the  trailing  edge,  and  to  be  aligned  with 
the  freestream  tnflnlty,  and  the  usual  boundary  conditions  of  tangential  flow 
apply  over  the  airfoil  and  sheet. 

Spence  solved  the  problem  for  the  specific  case  of  a flat  plate  airfoil 
and  simple  hinged  flap,  with  the  jet  injected  parallel  to  tf:u  flap.  These 
results  have  always  been  of  great  Interest  since  they  serve  as  a comparative 
baseline  for  assessing  the  utility  of  more  approximate  methods.  Spence's 
results  are  fairly  accurate  for  many  cases,  and  in  fact  give  good  results  well 
outside  the  normal  bounds  of  linearized  theory.  At  the  same  time,  titey  are 
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difficult  to  extend  to  include  real  wing  characteristics  such  as  a camber, 
large  leading  edge  devices,  multi-hinged  flap  segments,  and  momentum 
injection  which  is  nonplanar  with  the  flap. 

The  extension  of  such  concepts  to  three  dimensional  wings  immediately 
leads  to  far  greater  complications.  As  in  the  two  dimensional  case  the  jet 
sheet  vortex  loading,  with  its  interaction  with  the  wing,  is  the  critical 
driving  function,  so  that  the  impact  of  three  dimensionality  on  its  shape 
requires  careful  assessment.  The  jet  sheet  loading,  being  three  dimensional, 
also  sheds  its  spanwise  gradient  of  lateral  vorticity,  so  that  the  accumulated 
trailing  vorticity  increases  to  its  asymptotic  value  at  downstream  infinity. 


This  problem  has  been  attacked  in  numerous  ways.  Maskell  and  Spence(3) 

. developed  first  order  trends  by  collapsing  the  chordwise  dimension  of  both 

wing  and  jet  sheet  to  an  ellipticaily  loaded  lifting  line,  and  by  assuming  a 
flat  wake  with  the  usual  doubling  of  downwash  from  the  loaded  line  to  the  far 
field.  In  applying  Spence's  previously  developed  section  properties ( - ' an 
additional  factor  must  be  taken  into  account,  namely,  that  at  the  lifting  line 
the  downwash  is  aj,  but  that  the  jet  sheet  far  downstream  is  subjected  to  a 
downwash  angle  = 2<xj.  The  total  jet  turning  angle  is  therefore  (9 j - ctj)  - 
k (3OT-ctj),  where  0j  is  the  injection  angle  measured  from  the  freestream  direc- 

tion. The  jet  sheet  loading,  and  therefore  section  properties,  are  adjusted  to 
reflect  a curve  freastream.  The  section  properties,  in  local  flow  coordinates, 
are  therefore  in  part  determined  by  three  dimensional  downwash.  The  basic 
elliptic  loading  results  by  Maskell  and  Spence  can  give  good  lift  results  for 
nonelllptic  cases,  provided  that  appropriate  part  span  factors  are  used,  but 
they  are  not  suitable  for  calculating  span  load  distribution  and  induced  drag. 


Other  researchers  have  produced  workable  techniques  while  addressing 
k fairly  generalized  planforms  and  distributions  of  flap  chord,  flap  deflection, 

local  momentum  coefficients  and  injection  angle.  Some  of  the  more  notable 
i works  are  due  to  Das'**',  Lopez  and  Shen'5),  llssaman'^),  Hackett  and  lyman(7)t 

and  Davenport vo) . The  approach  reported  in  wh i s document  Is  in  many  ways  a 
lifting  line  counterpart  of  the  Hackett-lyman  approach,  but  with  some  modlfl- 
( cations  in  the  assumptions.  Davenport's  method,  which  became  available  to 

these  authors  fairly  late  in  the  study,  is  of  considerable  Interest  since  it 
is  also  a lifting  line  technique,  but  which  differs  somewhat  in  approach, 
assumptions,  and  method  of  solution. 

It  Is  pointed  out  that  the  sheet  vortex  loading  is  highest  in  the 
vicinity  of  the  wing,  followed  by  a rapid  decay.  Generally,  at  some  nominally 
"short"  distance  downstream,  the  substantially  one-to-one  conversion  of  jet 
| momentum  lift  to  sheet  vortex  lift  is  essentially  complete,  particularly 

Insofar  as  its  lift  interference  with  the  wing  is  concerned.  Accordingly,  the 
system  may  be  suitably  truncated  for  calculation  of  Interference  between  the 
jet  sheet  and  wing.  Moreover,  for  most  cases  of  present  Interest  the  sheer 
roll-up  occurring  between  the  wing  and  truncation  point  Is  Ignored,  although 
this  can  b«  significant  when  spanwise  gradients  are  sufficiently  large. 
Additionally,  questions  arise  concerning  the  nonplanar  character  of  the  wake, 
since  It  Is  outside  the  wing  reference  plane  and  it  may  have  considerable  warp 
* due  to  spanwise  variations  in  both  injection  angle  and  downwash  field.  Even 

so,  the  analysis  has  shown  that  a completely  planar,  linearized  approach  Is 
capable  of  giving  good  results  although  several  factors  neglected  in  the 
linearization  are  known  to  be  Individually  significant. 
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3.  APPROACH 


The  approach  taken  here  is  based  on  a completely  planar  system  using 
linearized  theory  and  linearized  boundary  conditions.  Moreover,  it  is  based 
on  the  assumption  that,  as  far  as  its  Influence  on  the  total  lift  of  the  wing 
is  concerned,  the  jet  sheet  can  be  approximated  by  a flat  plate  which 
possesses,  at  the  spanwise  station  considered,  the  inclination  of  the  injected 
jet  momentum  and  a (unknown)  length  in  the  streamwise  direction  which  is 
determined  by  properties  of  the  jet  sheet.  Details  will  be  discussed  later. 
The  solution  involves  iteration  of  the  spanwise-varying  length  of  a flat  sheet 
whose  inclination  is  that  of  the  Injected  momenta,  in  contrast  to  the 
classical  problem  of  shape  calculation. 

The  primary  objective  of  this  program  was  the  calculation  of  rolling 
moments  and  sidewash  in  sideslip.  However,  such  calculations  necessarily 
require  accuracy  In  span  loadings  and  a clear  distinction  between  vortex 
forces  and  "residual1'  efflux  momentum  forces,  thereby  providing  the  Informa- 
tion essentia)  to  an  accurate  assessment  of  Cl  - a and  Cl  ■ Cx  relationships. 
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SECTION  II 


THEORY 


The  analysis  procedure  given  in  this  report  consists  of  a swept  lifting 
line  approach  to  the  analysis  of  jet  flap  type  systems  in  sideslip.  The 
method  makes  use  of  the  Equivalent  Mechanical  Flap  Concept  (EMF)  for  calcula- 
tion of  powered  section  properties.  First  order  fuselage  interference  is 
taken  into  account.  There  are  certain  elements  and  special  considerations 
which  require  explanation.  The  following  text  is  confined  primarily  to  these 
areas  and  to  the  manner  in  which  they  fit  into  the  total  picture. 

This  section  considers  the  two  dimensional  EMF  in  both  uniform  and  curved 
streams,  the  particular  boundary  conditions  ,;ed,  fuselage  interference,  the 
nonlinear  system  of  equations  and  method  of  juration,  the  circulation  and 
residual  jet  lifts,  the  assessment  of  axial  force  {partly  as  a point  of  vali- 
dation for  span  load  distributions),  and  sidewash  and  downwash.  These  areas 
are  summarized  in  Figure  1. 


1.  THE  TWO  DIMENSIONAL  EQUIVALENT  MECHANICAL  FLAP  IN  A UNIFORM  FREESTREAM 
a.  Physical  Considerations. 

A typical  2-0  jet  flap  arrangement  is  illustrated  in  Figure  2.  The 
condition  for  dynamic  equilibrium  of  the  sheet  is  that  the  pressure  differen- 
tial across  the  jet  locally  balances  the  centrifugal  force  of  the  jet  efflux. 
For  small  perturbations  this  leads  directly  to 
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so  that  the  jet  sheet  lift  exerted  over  the  region  between  the  trailing  edge 
and  some  arbitrary  location  downstream,  say  x « x',  is 
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where  0 Is  the  total  injection  angle  relative  to  -he  freestream  direction. 
Thus,  the  accumulated  vortex  lift  equals  the  difference  in  the  lift  component 
of  the  jet  between  the  two  stations.  Typically  the  curvature  z"  is  large  near 
^ the  trailing  edge,  followed  by  a rapid  downstream  decay,  so  that'  Me  typically 

becomes  very  small  a short  distance  downstream.  This  is  especially  true  for 
nominal  momentum  coefficients.  The  distance  in  chords  required  for  §/9  to 
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degenerate  to  a specified  level  increases  with  increasing  momentum, 

Cp.  It  also  Increases  with  increasing  levels  of  the  power-off  lift  of  the 
mechanical  system,  since  a pre-existing  local  flow  inclination  of  the  trailing 
edge  reduces  the  effective  jet  disturbance  for  a specified  injection  angle. 
Thus  the  sheet  geometry  is  sensitive  to  Cp,  angle  of  attack,  flap  chord,  flap 
deflection,  airfoil  camber,  and  injection  angle. 

In  1958  Jacobs  and  Paterson^)  proposed  that  the  two  dimensional  analysis 
be  simplified  by  replacing  the  sheet  by  a finite-length  curved  hard  surface 
flap  characterized  by  semi-empirical  logarithmic-type  shape  functions,  and 
establishing  by  iteration  a chordwise  extent  of  flap  sufficient  to  carry  the 
total  vortex  lift  on  the  sheet.  The  results  of  their  study  compared  very  well 
with  experimental  data. 

Parallel  studies  considered  a simpler  case  in  which  the  sheet  was 
replaced  bv  a straight  mechanical  flap  which  was  deflected  parallel  to  the  jet 
injection  angle.  The  procedure  is  basically  simple,  since  it  is  only  required 
to  establish  with  thin  airfoil  theory  an  extended  length  Ac  which  carries  a 
lift  Cp  sin  6 in  the  presence  of  the  actual  airfoil.  The  calculated  lift 
results  were  also  quite  good.  This  is  the  approach  used  in  this  report. 
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b.  Derivation. 

For  completeness,  the  discussion  which  follows  reiterates  the 
origin*-  derivation  of  Jacobs  and  Paterson(9)  and  the  closed  form  solution 
given  uo  e recently  by  Hackett'?), 


U is  presumed  that  a flat  mechanical  flap  extension,  henceforth  referred 
to  as  the  Equivalent  Mechanical  Flap  (EMF),  can  be  used  to  represent  the  lift 
aspects  of  the  jet  sheet.  The  direct  jet  lift  and  accordingly  the  sheet 
vortex  lift,  is  Cp  sin  0,  and  an  EHF  length  will  be  established  so  that 
cu  sin  0 » 

In  order  to  do  this,  reference  is  made  to  Hunk's  integral  theorems  for 
thin  airfoils  (see  Ref.  10,  Section  E.11,  Equations  3.\k,  9-20).  These  are 
written 
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where  the  notation  is  adapted  to  that  shown  in  Figure  3.  These  equations 
combine  to  give 
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Consider  the  special  case  shown  in  Figure  3 where  dz/dx  is  constant  ove» 
an  element  of  airfoil  defined  by  0 < ip  < The  y distribution  due  to  that 
element  is,  based  on  the  preceding  expression. 
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The  Identity 
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allows  Eq  (I)  to  be  written 


) 

2V(-f) 


■/“ 


cot 


<P'  - ip  . 1 - cosip’ 

*-2  ♦ TsTSft-  dt 


which  is  directly  integrable  to 
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Now  consider  the  region  0 « * « v|,  which  wilt  subsequently  become  the  equi- 
valent mechanical  flap.  The  lift  induced  on  this  region  is 
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v;hich,  with  Eq  (2),  gives 
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By  rearranging  Eq  (3)  one  writes 
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which  Is  tb*  total  lift  In  the  range  0 < $ f vj . due  to  a constant  deflection 
-d*/dx  over  the  range  0 « $ * The  corresponding  total  airfoil  lift  Is 
obtained  from  Eg  (?)  with  • * (load  integrated  over  the  entire  chord). 
That  is,  with  -da/dx  s 4 
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= 2 (4 'i  + s I nij;)  6 

For  a section  lift  curve  slope  acg/aa  * 2it,  this  reduces  to  the  familiar  a6 
relationship 
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The  preceding  result  may  be  generalized  to  a Ac^  consisting  of  the  total 
contributions  of  N airfoil  components  (see  Hackett  and  Lyman  (7)).  However, 
the  present  case  is  restricted  to  a plane  airfoil  with  nose  flap,  a single 
main  flap,  and  the  EMF  as  show.’;  in  Figure 


The  elements  are  as  follows: 

Flap  Size 

Element  Parameter,  $ 

0 < tb  < it  it 


0 < |b  < tb 

- - n 

0 < ib  < ib 

- - O 

0 < <>  < 


*0 

*1 


Deflection  (-dz/dx) 
a - 6 

n 


6 

6 

6 


n 

f 
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The  elements  superimposed  to  give  the  airfoil  geometry  shown.  The 
following  geometric  relations  apply 
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The  length  of  the  equivalent  flap  (or  equivalent  ) Is  now  determined  by  the 
requirement  that  the  vortex  lift  exerted  on  It  equals  the  total  vortex  lift  of 
the  jet  sheet.  The  total  vortex  lift  on  the  EMF  is  obtained  by  applying 
Eq  (7)  to  all  wing  elements  listed  above  and  adding  the  results.  The  total 
vortex  lift  of  the  jet  sheet  (referred  to  chord  length  c')  is  given  by  the 
change  of  the  vertical  momentum  in  the  jet  sheet*  that  is,  by 


cy  Sin(ct  + 6f  + 6 j ) p- 


where  cy  is  referred  to  the  geometric  chord  c.  One  thus  obtains  the  final 
equivalence  relation 
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In  a similar  manner  the  total  airfoil  lift  is  obtained  by  superimposing 
the  contributions  of  the  four  elements  (sec  Eqs  (8),  (9)).  The  resulting  lift 
coefficient  based  on  physical  chord  c,  is  then 


\ m U T 7twx  [(°  “ *n)  * °«  (*n)  + Sf  * V**1  *j] 


(15) 


11 


Since  the  jet  is  ultimately  turned  back  to  the  direction  of  the  free- 
stream,  and  the  aerodynamic  drag  of  the  combined  wing  and  jet  sheet  is  zero, 
the  axial  force  is  given  by 


CX  = Cu  (16) 

as  stated  by  the  classical  thrust  recovery  arguments. 

Calculation  of  <J»i  in  (14)  to  match  the  value  of  H specified  in  Eq  (13) 
required  use  of  the  relations  (11)  and  (12).  The  quantity  H is  single  valued 
in  ibj  over  the  range  0 < l^l  < tr(dH/dtj)i  >0),  so  that  the  initial  guess  i|ji  =ir/2 
establishes,  by  comparing  left  and  right  sides,  which  quadrant  ^ falls  within. 
The  second  guess  is  the  mid-range  of  the  appropriate  quadrant  (^i  = ttA  or 
<|>l  = 3“A)  > so  that  the  second  sector  containing  is  reduced  in  size  to  ir/2n 
in  n steps. 

The  iterated  value  of  ^ , with  the  corresponding  ip  and  i|>n,  are  used  in 
(15)  for  calculation  of  the  lift. 

It  is  reiterated  that  i|)lf  while  sensitive  to  the  direct  jet  lift  C£  sin 
(a  + 6f  + 6 j ) , is  also  sensitive  to  the  size  and  deflection  of  each  airfoil 
element.  This  is  of  course  reflected  as  a sensitivity  of  jet  induced  super- 
circulation lift  to  the  geometry  of  the  mechanical  system. 


2.  EFFECT  OF  LONGITUDINAL  VARIATION  IN  THE  INDUCED  ANGLE 

A longitudinal  variation  of  the  downwash  angle,  relative  to  its  value  a[ 
at  the  wing,  will  change  the  total  sheet  vortex  load  and  its  longitudinal 
distribution.  This  will  be  taken  into  account  only  to  the  extent  that  it 
changes  the  loading  induced  by  the  jet  on  the  physical  wing.  Only  the  forward 
portion  of  the  jet  sneet  need  be  taken  into  account  for  determination  of  the 
lift.  The  portion  of  the  sheet  taken  into  account  is  subjected  to  the  induced 
angle  aj  at  its  upstream  edge  and  to  the  induced  angle  3 at  its  downstream 
edge  (see  Figure  5)-  The  change  of  jet  momentum  in  the  direction  normal  to 
the  local  flow  at  the  wing  is 

cy  [sin (a  - aj  + 6f  + 6j)  - sin  (3  - ctj)] 

which  is  ttn  revised  vortex  lift  to  be  carried  by  the  EMF. 

The  appropriate  lift  equations  are 

jj-  cy  [sin (a  - a;  + + 5j)  - sin(3  - aj)]  = H (a  - aj,  i|»i)  (17) 


where  H(a  - aj,  t^)  is  obtained  from  Eq  (15)  with  a replaced  by  a - aj. 
Correspondingly,  section  circulation  lift  is  obtained  from  Eq  (15)* 
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CZ  ~ c£v  cosaj  + Cp  sin  S 


(19) 


~ C£y  + cu  sin  8 


It  is  emphasized  that  8 - a;  is  the  variation  in  downwash  angle  across 
that  portion  of  the  jet  sheet  which  is  presumed  to  determine  the  super- 
circulation or  jet  interference  lift  on  the  wing. 

It  is  pointed  out  that  8 = a;  at  the  wing.  If  the  entire  sheet  is  taken 
into  account,  then  8 = where,  for  an  elliptically  loaded  system,  8 = 2aj 
(flat  wake)  and  8 ==  (4/7t2)ct ? (rolled  wake).  Thus,  8 - a]  is  at  most  of  order 
aj  (flat  wake)  or  -(l  - k/n2)  a;  (rolled  wake). 

Since  the  variation  8 - aj  is  proportional  to  aj  it  is  significant  only 
when  a;  itself  is  large.  This  condition  is  best  typified  by  large  Cu  and 
concentrated  part  span  blowing  (high  lift,  low  effective  aspect  ratio).  For 
such  extreme  cases  consistently  good  lift  results  have  been  obtained  simply 
by  assuming  that  8 a a;  and  neglecting  the  downwash  variation  entirely.  In 
the  computer  programs  included  in  this  report  8 is  set  equal  to  a}. 


3.  THE  LIFTING  LINE  METHOD  FOR  SPAN  LOAD  CALCULATION 

a.  The  Boundary  Condition 

The  method  reported  here  uses  the  strip  approach  and  simple  sweep 
theory.  Most  methods  based  on  this  approach  Impose  the  boundary  condition  of 
tangential  flow  at  the  three-quarter  chord  location.  A strict  interpretation 
required  the  section  lift  curve  slopes  to  be  2u;  otherwise,  a chordwise 
adjustment  in  control  point  location  Is  required.  In  the  present  case  the 
lift  curve  slope  of  the  airfoil  plus  EMF  Is  taken  to  be  2rr,  based  on  the 
effective  chord  c'.  However,  the  effective  chord  length  c*  is  itself 
initially  unknown  in  the  3"0  field  since  It  is  dependent  on  the  Induced  angle, 
aj.  Thus,  enforcement  of  boundary  conditions  at  the  3/A  chord  point  requires 
that  the  physical  locations  of  all  control  points  be  simultaneously  iterated 
along  with  ot|.  This  would  further  require  that  Influence  coefficients  be 
recalculated  in  each  iteration  step.  (See  Figure  6a.) 

In  order  to  minimize  both  computer  time  and  potential  convergence 
problems  the  following  approach  is  used.  It  is  assumed  that  section  proper- 
ties are  merely  the  response  to  an  effective  angle  of  attack  ae  « a - aj, 
where  a}  is  evaluated  at  the  lifting  line  Itself,  normally  taken  as  the  1/A 
chord  line.  The  boundary  condition  to  be  met  at  specified  control  points  is 
then  ae  * a - aj.  The  vortex  or  circulation  lift  due  to  the  combined 
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mechanical  airfoil  and  EMF  (c for  the  equivalent  airfoil)  is  represented  as 
a bound  vortex  on  the  quarter  chord.  The  gradient  of  the  distribution  of 
bound  vorticity  establishes  the  strength  of  the  trailing  vortices  in  the  usual 
manner.  The  induced  angle  consists  of  that  induced  by  the  trailing  vortices, 
and,  when  the  wing  is  swept,  by  the  contribution  of  the  bound  vortex  repre- 
senting the  opposite  wing  panel.  (See  Figure  6b.)  This  is  in  contrast  to  the 
3/4  chord  boundary  condition  which  requires  the  downwash  contribution  from  all 
vortex  elements. 

The  computational  procedure  developed  in  this  study  is  based  on  a straight 
lifting  line  from  wing  root  to  tip,  onto  which  all  section  properties  are 
lumped.  In  practice  the  location  of  a constant  percent  chord  line  is  shifted 
due  to  either  mechanical  flap  extension  or  the  EMF  extension.  This  shift 
depends  upon  the  spanwise  location.  For  a constant  longitudinal  shift  of  the 
lifting  line,  all  equations  remain  the  same  since  all  properties  are  lumped 
onto  the  line  itself,  in  this  sense  the  procedure  has  a limited  self-adaptive 
capability.  The  user  may  elect  to  use  any  constant  sweep  line  which  is 
believed  to  most  effectively  represent  the  wing  and  EMF  by  a straight  bound 
vortex.  Good  results  have  been  obtained  by  using  the  quarter  chord  of  the 
basic  wing  planform,  which  indicates  that  this  selection  is  not  critical.  o 


b.  Influence  Coefficients  of  Lifting  System 

The  loaded  line  and  its  associated  trailing  vortex  system  is 
identical  to  the  Weisslnger  approach  given  in  Ref.  11  with  the  exception  that 
influence  coefficients  for  downwash  are  calculated  at  the  lifting  line  itself 
rather  than  at  the  3/4  chord  location. 

The  calculation  of  downwash  influence  coefficients  is  a straightforward 
modification  of  those  In  Ref.  (11),  and  only  summary  comments  follow.  The 
method  considers  a distribution  of  bound  vortex  strength,  r which  Is  continu- 
ous along  the  lifting  line  and  which  vanishes  at  each  wing  tip.  The  T 
distribution  is  allowed  to  be  unsymmetrical . The  downwash  integrals  include 
both  T (bound  vortex)  and  dT/dy  (trailing  sheet)  in  the  integrand.  Integration 
by  parts  allows  the  complete  downwash  Integral  to  be  written  as  a function  of 
geometry  and  dr/dy  only.  By  writing  r In  terms  of  p (y  ■ b/2  cos  b)  and  using 
Multhopp's  quadrature  formula,  the  final  Influence  coefficients  aVfl  are 
obtained,  giving 


m 

°lv  ■ *vn  Gr«»  v - 1,  2 

n*1 


m 


(20) 


where  afv  Is  the  angle  Induced  in  the  streamwlsc  plane  at  span  station  and 
Gn  ■ r(pn)/bV  represents  the  bound  vortex  strength  at  span  station  pn. 

Stations  and  are  standard  for  each  specified  m number  of  control  points, 
that  is,  $v  * v»/(m+l). 


For  a given  geometry  the  influence  coefficients  avn  can  be  computed  once 
and  for  ail.  They  do  not  depend  upon  the  sizing  of  the  EMF,  which  changes  in 
each  iteration  step.  The  iteration  is  described  later. 


c.  Fuselage  interference  on  the  Wing 

The  first  order  effect  of  the  fuselage  on  the  wing  is  generation  of 
a local  flow  inclination  due  to  both  pitch  and  sideslip,  which  is  treated  here 
by  superimposing  onto  the  wing  an  equivalent  change  in  local  geometric  angle 
of  attack.  The  method  used  is  extremely  simple,  in  that  the  wing  is  assumed 
to  be  submerged  in  the  field  generated  by  the  isolated  body.  Additionally, 
the  body  is  represented  by  a cylinder  which  extends  to  infinity  in  both 
directions,  and  whose  constant  cross  section  is  chosen  as  that  of  the  body  in 
the  vicinity  of  the  inboard  wing,  where  the  effects  are  significant.  For 
simplicity,  the  body  cross  section  is  approximated  by  an  ellipse,  of  arbitrary 
axis  ratio,  whose  center  is  arbitrarily  located  in  the  vertical  direction.  The 
calculation  of  the  fuselage  cross  flow  is  nonlifting  about  the  fuselage  itself, 
so  that  its  impact  is  confined  to  wing  interference.  Therefore,  in  the  calcu- 
lation program  the  interference  angle  of  attack  on  the  wing  is  treated  as  an 
equivalent  geometric  wash-in  or  wash-out  which  is  constant  during  the  span 
load  iteration.  In  all  subsequent  discussion  the  geometric  angle  of  attack  of 
each  airfoil  section  is  understood  to  consist  of 
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where  the  last  term  is  the  result  of  geometric  dihedral  in  sideslip.  The 
calculation  of  Aa/agocjy  and  Ao/8  is  given  in  Appendix  A. 


d.  The  Nonlinear  System  and  Method  of  Iteration 

t 

(1)  Simple  Sweep  Considerations.  In  keeping  with  simple  sweep 
theory  the  boundary  condition  (aeaci-'*|)  is  Imposed  in  the  simple  sweep  plane. 
However,  the  equation  is  unchanged  in  streamwlse  coordinates  since  the  trans- 
formation Is  the  same  for  all  a's.  The  system  of  equations  in  streamwlse 
coordinates  is  therefore 
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Equation  (22)  represents  a system  of  m equations  for  m unknown  values  of  aev» 
since  av  and  aVf1  are  specified  by  the  geometry  and  since  Gn  is  ultimately 
expressible  in  terms  of  aen.  The  circulation  function  Gn(aen)  ,s  expressible 
in  terms  of  the  section  circulation  lift  coefficient  given  by  Eq  (18),  with 
a - a;  replaced  by  ae  (note  that  sweep  must  be  taken  into  account).  The 
quantity  ipj  which  occurs  in  Eq  (18)  is  expressed  by  Eq  (17)  in  terms  of 
geometric  quantities,  momentum  coefficient,  and  ae.  Thus,  one  can  write 


Gn  = Gn(aen) 


(24) 


where  this  relation  encompasses  a rather  lengthy  algorithm  involving  Eqs  (18) 
and  (17),  and  the  simple  sweep  conversions  given  in  Figures  7 and  8. 

(2)  Iteration.  For  the  special  case  of  zero  momentum  coefficient 
the  quantity  ipj  is  zero  so  that  Eq  (22)  is  linear  in  the  unknowns  aev-  In 
the  general  case  the  system  is  nonlinear  and  iteration  is  required. 

Equation  (22)  may  be  written 
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or,  in  vector  form 


(26) 


where  the  star  superscript  indicates  the  actual  solution. 

In  the  vicinity  of  the  solution  ae  ■ ae*  , cy  ■ 0,  one  may  make  the 
linear  approximation 
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or  alternatively 
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Now  m3ke  a guess  for  the  m values  of  aey.  These  are  used  In  (25)  to  calculate 
the  error  ev.  The  root  mean  square  error  is  calculated  for  future  reference. 
The  guessed  ae  and  the  calculated  error  e are  substituted  into  Eq  (28).  The 
result  Is  m number  of  equations  which  are  linear  in  the  m values  of  ae*  and 
linear  in  the  m2  number  of  partial  derivatives  9aev/9en*  The  r.i  equations 
involve  a total  of  m (m+1)  unknowns.  By  making  a total  of  m + 1 linearly 
independent  guesses  one  has  m (m  + 1)  equations  and  m (m  + 1)  unknowns. 

The  equations  are  solved  for  ae*  which  are  then  used  In  (25)  to  obtain 
corresponding  errors  and  the  RMS  error.  The  guess  whose  RMS  error  is  greatest 
is  replaced  by  the  "improved  guess"  ae**  The  approximately  linear  relation- 
shlp(28)  is  thus  updated  and  a new  value  of  ae*  is  calculated.  The  Iteration 
is  terminated  when  the  RMS  value  of  ae*  Is  within  a prescribed  tolerance. 


e.  Three  Dimensional  Residual  and  Total  Lift 

With  the  iterated  effective  angle  of  attack  known  at  m points  across 
the  span,  the  corresponding  section  circulation  lift  is  calculated  and  numeri- 
cally Integrated  to  give  total  circulation  lift  for  the  wing  and  EMF  combined. 
There  is,  in  addition,  a residual  momentum  lift,  which  must  be  taken  into 
account.  Within  the  framework  of  slmpte  sweep  theory  (infinite  yawed  wing), 
there  are  no  pressure  or  area  gradients  in  the  spanwise  direction,  so  that 
loads  and  momentum  changes  are  due  only  to  the  velocities  and  geometry  which 
project  onto  the  simple  sweep  plane.  Although  particles  of  fluid  move  along  a 
nonplanar  path,  the  pressure  distribution  is  defined  by  the  propagation  of 
that  motion  onto  the  simple  3weep  plane.  The  classical  arguments  are  readily 
extended  to  include  the  momentum  injection  when  the  momentum  vector  lies  in 
the  simple  sweep  plane. 

In  the  general  case  the  jet  has  a component  of  momentum  normal  to  the 
simple  sweep  plane,  and  the  skew  angle  between  the  jet  and  the  simple  sweep 
plane  is  denoted  by  e as  shown  In  Figure  8.  The  component  of  momentum 
Jn  " J which  lie*;  in  the  simple  sweep  plane.  Is  turned  back  parallel  to 

the  projection  of  the  Incident  stream  onto  the  simple  sweep  plane.  The  compo- 
nent which  Is  Injected  perpendicular  to  the  simple  sweep  plane,  jA  » j sine, 
is  unaffected  and  It  must  be  transported  Intact  during  the  turning  process. 
Thus,  the  simple  sweep  theory  still  applies  to  this  case  without  modification. 
Therefore,  the  projections  of  all  relevant  quantities,  including  the  injected 
momentum  vector,  onto  the  simple  sweep  plane  have  been  used  as  the  basis  for 
calculation  of  section  circulation  properties. 
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Consider  now  the  corresponding  residual  lift.  This  lift  teim  consists  of 
the  projection  of  two  vectors  onto  the  airplane  lift  axis.  The  two  vectors 
are  (a)  the  component  which  is  injected  normal  to  the  simple  sweep  plane  and 
transported  intact  during  the  turning  in  the  simple  sweep  plane  and  (b)  the 
airplane  component  which  has  been  turned  parallel  to  the  incident  stream  in 
the  simple  sweep  plane.  The  geometry  Is  shown  in  Figure  9 The  components 
perpendicular  to,  and  contained  by  the  simple  sweep  plane  are  denoted  by  j^ 
and  jn,  respectively.  The  residual  lift  is  then 


^residual  = cos(A  1 + k s]n('h  1 sin  ais  ^ 

where  the  upper  and  lower  signs  correspond  to  the  left  and  right  hand  panel:., 
respectively.  In  terms  of  the  injection  skew  angle,  e,  one  has 

jn  = j cose 

jA  - j sine 


from  which  the  residual  lift  , Eq  (29),  becomes 

Vesidual  = j cos  (A  ± B - e)  sin  ajs  (30) 

This  expression  is  integrated  across  the  span  to  give  the  total  residual  lift. 

It  is  noted  that  the  residual  lift  in  Eq  (30)  is  not  defined  exclusively 
by  conditions  In  the  simple  sweep  plane.  At  6f  + 6j  a 0 the  preceding  term 
corresponds  exactly  to  the  momentum  component  which  is  injected  onto  the  plane 
containing  the  freest  ream  velocity  vector  and  the  airplane  lift  axis  (pitch 
plane  when  B * 0) . In  general,  for  5f  + 6j  ¥ 0,  the  component  injected  within 
the  aforementioned  plane  is 

jp  » J(cos2(A  t B -e)  + sin2(6f  + 6j)  (cosze  sin2(A  + S)) 
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-2  [(1  - cos(6f  + 5j))  sine  cose  sln(AtB)  cos(A  i 8)  r 
If  6f  + 6j  is  not  too  large,  this  may  be  approximated  by 
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In  the  ranges  of  A + B and  6f  ♦ for  which  linearized  theory  is  expected  to 
apply,  this  expression  is  not  seriously  different  from  the  term  j cos(A$8-c) 
used  in  the  calculation  of  the  residua)  lift. 


18 


f.  Rolling  Moments  In  Sideslip 


Rolling  irtoments  in  sideslip  are  obtained  by  Integrating  the  moment 
of  the  total  lift  distribution,  regardless  of  how  direct  jet  thrust  lift  is 
split  between  residual  and  vortex  lift.  This  is  given  by 

1 
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The  present  computer  program  is  restricted  to  symmetrical  geometry  and 
momentum  injection  at  arbitrary  sideslip.  The  extension  to  cover  the  more 
general  case  of  engine-failed  rolling  moments  and  aileron/flap  lateral  trim 
and  control  is  straightforward  since  only  a c^  and  aileron/flap  input 
generalization  is  required. 

g.  Axial  Force 

This  portion  of  the  study  was  originally  conducted  as  a check  on  the 
total  span  load  distribution,  its  split  between  residual  and  vortex  lift,  and 
the  entire  concept.  The  program  modifications  to  do  this  are  negligible.  The 
Cx  calculations  to  date  give  additional  confirmation  of  both  concept  and  com- 
putational accuracy  and  are  of  course  extremely  useful  in  themselves. 

The  following  considerations  are  pertinent.  First,  lift  is  of  first 
order  in  angle  of  attack,  whereas  axial  force  is  of  second  order.  Second, 
lift  is  obtained  as  the  sum  of  circulation  lift  and  jet  efflux  momentum  lift, 
whereas  axial  force  is  obtained  as  a difference  between  jet  efflux  momentum 
(thrust)  and  a circulation  (Induced)  drag. 

Inasmuch  as  a redirection  of  the  jet  produces  changes  In  the  total  level 
of  circulation,  it  follows  that  axial  force  Involves  differences  in  terms 
which  can  individually  involve  terms  of  second  order  Ir.  angle  of  attack,  where 
such  second  order  thrust  and  drag  terms  become  Increasingly  interactive  for 
increasing  levels  of  aeropropulslve  integration. 

Specifically,  it  has  been  assumed  in  earlier  discussion  that  the 
character  of  the  remote  wake  has  negligible  influence  on  the  overall  lift 
characteristics.  It  is  also  assumed  that  the  local  jet  Is  bent  back  parallel 
to  the  local  downwash  angle  and  that  additional  jet  bending  further  downstream 
is  sufficiently  remote  that  its  effects  on  lift  can  be  Ignored.  Such 
additional  turning  results  in  an  additional  conversion  of  thrust  lift  into 
sheet  circulation  lift.  This  conversion  is  on  a one-to-one  basis  since  it  is 
assumed  to  produce  no  significant  effect  on  the  wing.  The  corresponding 
changes  in  residual  thrust  and  in  circulation  or  vortex  drag  do  not  compensate 
because  of  their  second  order  nature. 

The  axial  force  coefficient  is  calculated  from  Trefftz  plane  considera- 
tions, for  which  a more  precise  breakdown  of  total  circulation  and  thrust 
forces  is  required.  The  bookkeeping  difference  is  Illustrated  in  Figures  10(a) 
and  10(b).  The  wing  circulation  lift,  labeled  "A,"  is  identical  In  both  the 
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near  and  remote  fields.  The  total  jet  efflux  lift,  consisting  of  a circula- 
tion lift  "B"  and  a residual  lift  "C",  is  also  identical.  The  particular 
breakdown  of  "B"  and  "C"  corresponding  to  the  near  field  calculation  (Figure 
10(a)  is  not  proper  for  axial  force  calculation.  The  more  correct  far  field 
breakdown  is  shown  in  Figure  10(b). 

The  method  outlined  In  Appendix  B takes  into  account  this  finer  break- 
down while  making  use  of  shape  factors  based  on  the  spanwise  distribution  of 
circulation  forces  obtained  from  the  near  field  lift  calculation  already 
described. 


h.  S! dewash  and  Oownwash 

In  sideslip  the  planview  of  the  trailing  vortex  system  is  skewed 
from  the  airplane  center, ine  by  the  sideslip  angle  B (approximately).  This 
lateral  displacement  Is  of  course  taken  into  account  in  sidewash  and  downwash 
calculation.  However,  the  effects  of  the  wake  skew  on  the  influence 
coefficients  avn  is  small  and  therefore  neglected  in  the  lift  calculation.  The 
lift  calculation  establishes  the  vortex  strengths  from  which  sidewash  and 
downwash  are  calculated  according  to  the  Biot-Savart  law. 


k.  SUMMARY  OF  METHOD 

The  preceding  discussion  has  dealt  with  the  individual  program  elements, 
special  considerations,  and  how  they  are  used  in  the  overall  analysis. 

in  order  to  give  a clearer  overall  view,  the  program  with  its  more 
essential  features  is  summarized  In  the  functional  flow  diagram  shown  In 
Figure  11.  The  input-output  details,  with  an  example  problem  and  print-outs, 
are  given  in  Appendix  0.  Computational  details  are  not  shown,  but  can  be 
obtained  if  required  from  the  computer  program  listing  given  in  Appendix  E. 
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SECTION  III 


COMPARISON  OF  THEORY  AND  EXPERIMENT 


The  comparisons  shown  In  this  section  illustrate  the  capabilities  of  the 
method.  They  are  typical,  are  in  no  way  selective,  and  are  generally  quite 
good  where  the  theory  is  expected  to  apply.  Poor  lift  comparison  is  generally 
confined  to  normal  upper  surface  stall,  to  lower  surface  separation  on  leading 
edge  devices  at  negative  angles  of  attack,  and  to  cases  where  Cy  is  too  small 
to  maintain  flow  attachment  over  the  flaps.  Rolling  moment  comparisons  are 
good  except  for  cases  involving  large  pylons  and  nacelles,  which  are  not  con- 
sidered in  this  analysis 


t.  EASIC  LIFT  COMPARISONS 

Most  comparisons  shown  here  were  made  early  In  the  program  development 
when  unsymmetrical  cases  and  axial  forces  were  not  yet  treated,  and  are  shown 
mainly  to  illustrate  the  overall  accuracy  of  the  method. 

Figure  12  shows  a basic  lift  comparison  for  a rectangular,  AR  = 2.75  wing 
without  a mechanical  flap.  The  data  are  from  Reference  12  (Williams  & 
Alexander).  It  is  noted  that  Das  tested  a substantially  identical  model 
(Reference  13)  for  which  the  comparison  is  much  better,  as  shown  in  Figure  13- 

The  lift  curve  slope  at  a = 0 is  shown  as  a function  of  Cu  in  Figure  Ik 
for  an  AR  » 5 rectangular  wingOM.  The  three  sets  of  data  shown  correspond 
to  different  freestream  velocities,  for  which  the  differences  are  minor.  The 
theoretical  curve  does  not  consider  this  since  it  is  a function  of  Cu  only. 

The  lift  at  a • 0 is  shown  in  Figure  15,  as  a function  of  Cu,  for  both  In- 
board and  outboard  half  span  blowing.  The  theory  correctly  shows  the 
Increased  effectiveness  of  inboard  blowing. 

Figure  16  shows  a comparison  for  a more  realistic  I BF  system  tested  by 
BoelngOS).  The  blown  flap  extends  outboard  to  n ■ 0.75  and  has  a AO  degree 
deflection.  An  additional  Jet  deflection  angle  of  about  10°  occurs  because  of 
the  upper  surface  trailing  edge  angle.  The  comparison  is  extremely  good. 

Figure  17  Illustrates  the  major  features  of  a large  scale  I BF  model 
tested  In  the  NASA-Ames  kO  x 80  foot  wind  tunnel  as  a Joint  AFFDl/NASA-Ames/ 
Lockheed  program(l6).  The  major  features  are  segmented  aileron-flap  blowlhg, 
the  Jacobs-Hurkamp  expanding  duct  arrangement,  and  the  trailing  edge  control 
flap.  The  control  flap  has  the  capability  of  modulating  the  Jet  injection 
angle  at  a fixed  flap  setting.  Some  lift  comparisons  are  shown  in  Figures  18 
through  23  as  various  flap  deflections  and  blowing  arrangements.  Lift  perform- 
ance with  full  span  30®  flaps  and  full  span  blowing  (constant  section  Cu  is 
shown  in  Figure  18.  The  effect  of  increa.ing  inboard  flaps  to  60  degrees  is 
seen  by  comparing  with  Figure  19.  The  60/30  flap  arrangement  of  Figure  19  is 
shown  In  Figure  20,  for  which  Jet  flap  blowing  is  restricted  to  the  60  degree 
inboard  flap  with  8LC  level/  only  on  the  aileron.  For  these  cases  the 
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comparison  is  generally  good.  Figures  21  and  22  show  Cl  and  Cl  at  a = 0 as 
functions  of  Cy  for  full  span  blowing  with  inboard/outboard  deflections  of 
60/30  degrees  and  60/60  degrees,  respectively.  It  is  clearly  seen  that 
excellent  comparison  is  reached  at  high  blowing  levels  typical  of  the  design 
condition.  It  is  noted  that  the  blowing  split  between  flap  knee  and 
trailing  edge  is  13  vs.  87  percent  of  the  total  flap  Cy. 

Figure  23  compares  theory  and  experiment  for  inboard/outboard  flaps  at 
60/30  degrees,  with  inboard  jet  flap  blowing,  and  with  the  control  jet  de- 
flected 40  degrees.  The  jet  Is  thus  Injected  at  100  degrees  relative  to  the 
wing  plane,  so  that  the  optimistic  theoretical  result  is  to  be  expected. 


2.  LIFT  AND  AXIAL  FORCE  FOR  ANALOGOUS  JET  FLAP  SYSTEMS 

This  section  compares  calculated  lift  and  axial  force  with  experiment  for 
four  large  scale  models  tested  In  the  NASA-Ames  40  x 80  foot  tunnel.  The 
systems  are  the  Internally  Blown  Flap  discussed  previously,  the  Augmentor  Wing, 
Externally  Blown  Flap  and  Upper  Surface  Blowing  (References  16,  17,  1 8 , 19). 

The  basic  clean  planforms  and  overall  arrangements  are  similar  for  the  four 
ccses.  System  details  are  not  shown;  schematics  are  illustrated  In  Figures  17,  24 
2j  and  26.  The  EBF  and  USB  calculations  are  based  on  estimated  jet  spreading. 
Augmentor  wing  calculations  use  $CU,  rather  than  Cy,  where  <J>  is  the  augmenta- 
tion ratio  under  static  conditions  and  Cy  Is  based  on  primary  (ejector) 
thrust..  The  calculated  augmentor  wing  C*  values  shown  have  not  been  corrected 
for  flap  augmentor  ram  drag. 

The  compared  results  are  summarized  in  Figures  27  through  30  for  each  of 
the  four  high  lift  concepts  at  typical  takeoff  and  landing  configurations,  it 
is  clearly  seen  that  the  present  method  accurately  predicts  the  first  order 
perfo  -mance  differences  among  these  concepts.  These  calculations  shown  In 
Figures  27  through  30  are  the  first  obtained  by  this  method  by  persons  other 
than  the  authors.  They  were  obtained  by  Messrs.  Y,  T.  Chin  and  Lee  Brandt,  of 
the  Lockl.ced-Georgla  Company,  without  prior  experience  In  either  program  usage 
or  "optlmiiud"  geometric  interpretation  of  certain  input  parameters  (e.g. 
optimum  representation  of  USB  Coanda  flaps  by  a single  flop  element). 

A typical  span  loading  is  shown  in  Figure  Jt  although  no  experimental 
data  are  available  for  comparison.  This  corresponds  to  the  US3  configuration 
cited  in  the  preceding  section  at  Cy  • 1.2,  e ■ 5°,  and  with  flaps  in  landing 
position  (see  Figure  3<).  The  wing  circulation  lift,  EMF  circulation  lift 
and  residual  momentum  lift  are  indicated  in  the  figure.  The  load  gradients 
between  successive  flap  p&.iels  are  apparent,  it  is  especially  interesting 
that  the  fraction  of  direct  lift  converted  to  EHF  vortex  lift  i»  about  25 
percent.  This  is  the  direct  result  of  the  severe  load  gradients  at  the  edges 
of  the  blown  section,  which  in  this  case  cause  local  downwash  angles  of  30  to 
3$  degrees.  The  effective  injection  angle  (relative  to  the  local  flow)  is 
thus  reduced  by  about  half. 
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3. 


ROLLING  MOMENT  DUE  TO  SIDESLIP 


The  majority  of  comparisons  shown  here  are  for  the  large  scale  IBF  model 
shown  earlier  in  Figure  26(20).  These  tests  were  made  tail -on.  Calcu- 
lated wing-body  results  have  therefore  been  corrected  for  rolling  moments  due 
to  the  rather  large  vertical  tail.  Sidewash  calculations  were  unavailable 
vdien  the  comparisons  wer's  made,  so  that  the  vertical  tall  correction  was  based 
on  standard  DATCOM  technlques(2l) . This  In  no  way  invalidates  the  comparisons 
since  the  expected  differences  in  vertical  tall  rolling  moment  are  small  com- 
r pared  to  the  configuration  total.  Table  1 summarizes  the  range  of  variables 

t covered  in  the  sideslip  runs.  The  variables  covered  are  angle  of  attack, 

momentum  coefficient,  main  flap  deflection,  and  control  flap  deflection.  All 
runs  were  at  30  degrees  aileron  deflection,  with  jet  flap  levels  of  blowing 
over  the  main  flap  and  BLC  blowing  on  the  aileron. 

Rolling  moments  fc.'  uniform  30  degree  flap/aileron  deflection  are  shown 
in  Figures  32(a)  through  (c) , for  a » 5 and  12  degrees  and  for  Cu  * 0.51,  1.01 
and  1.52.  The  sensitivity  of  wing-body  rolling  moment  to  angle  of  attack  is 
compensated  by  the  shift  in  vertical  tail  contribut’on.  It  is  surprising  that 
the  calculations  are  rather  insensitive  to  variations  in  Cu;  nevertheless,  this 
fact  Is  supported  by  experimental  data. 

The  effect  of  increasing  the  main  flap  deflection  to  60  degrees  is  seen  by 
comparing  Figures  32(b)  and  32(d).  At  both  o * A and  12  degrees  this  effect 
increases  the  rolling  moment  by  about  35  percent. 

Figures  32(e)  through  (g)  cover  the  same  a and  Cj,  range  as  Figures  32(a) 
through  (c) . The  main  flap/a  I leron/control  jet  deflections  are  60/30/20/degrees, 
respectively.  A general  increase  in  rolling  moment  is  apparent  relative  to  the 
corresponding  30/30/0  degrees  arrangement. 

' Figure  32(h)  shows  results  at  ■ 1.05  only  for  the  &0/30/A0  arrangement. 

Comparisons  of  Figures  32(d),  (f)  and  (h)  show  a sizeable  increase  in  rolling 
moment  due  to  increasing  control  flap  deflection. 

/ 

Rolling  moment  is  shown  In  Figure  33  for  the  large  scale  USB  model  noted 
eariler(19).  The  wing-body  calculation  corrected  for  vertical  tall  compares 
less  favorably  than  did  comparable  IBF  results,  which  is  most  likely  due  to 
the  large  over- the -wing  nacelles.  A coarse  assessment  of  the  nacelle  contri- 
bution Is  made  as  follows.  The  effect  of  fuselage  size  scales  approximately 
as  (dlameter)2  provides  that  the  wing  and  body  retain  the  same  relative 
position;  in  this  case  the  fuselage  and  nacelles  are  both  approximately 
snugged  to  the  wing.  Thus, 


23 


for  which  It  has  been  assumed  that  the  fuselage  effect  and  nacelle  effects  are 
pure  couples.  For  Dnac/Dfus  * 0.7,  the  nacelles  are  expected  to  cancel  the 
fuselage  effect.  The  correction,  &Cj.g  - +0.001  per  degree,  gives  considerable 
improvement,  as  seen  in  Figure  33. 

The  last  comparison  is  given  in  Figure  3**  for  the  Boeing  I BF  model  noted 
previously  (see  Figure  16,  Reference  15).  The  poor  comparison  is  probably  due 
primarily  to  p/lon-nacel le  crossflow  blockage  in  sidewash.  The  assumption  that 
fuselage  crossflow  is  effectively  blocked  allows  removal  of  about  ACj^  * - 0.001 
per  degree.  The  cascade  or  flow  straightening  effect  of  the  nacellesPand 
pylons  also  reduces  the  effective  wing  sideslip.  Moreover,  the  pylon  sideloads 
due  to  sideslip  shed  vortices  at  each  wing  pylon  juncture.  The  resulting  four 
couples  give  an  additional  positive  increment  in  C&.  No  attempt  is  made  here 
to  quantify  the  nacelle-pylon  interference,  which  is  an  area  requiring  further 
work. 


SECTION  IV 
EFFECTS  STUDIES 


This  section  presents  a limited  number  of  effects  studies.  The  dependent 
variables  are  rolling  moment  and  sidewash  derivatives;  independent  variables 
are  fuselagt  shape,  size,  and  vertical  location  relative  to  the  wing,  wing 
sweep,  goemetric  dihedral,  some  planfortn  effects,  and  distribution  of  trailing 
edge  blowing  and  injection  angle. 

It  is  not  intended  to  provide  a wide  range  of  working  curves,  but  rather 
to  illustrate  the  utility  of  the  method  and  some  first  order  trends.  For  this 
reason  the  effects  are  generated,  in  general,  as  perturbations  to  a baseline 
configuration,  which  is  chosen  as  the  large  scale  I6F  model  discussed 
earl ier(20) . The  total  blowing  coefficient  (based  on  wing  reference  area)  is 
Cu  s 1.01,  consisting  of  local  blowing  levels  (based  on  extended  chord)  of 
1 .0294  for  sections  out  to  70%  span  and  0.2628  for  sections  from  70%  span  to 
the  wing  tip.  The  wing  has  a full  span  leading  edge  device  which  contributes 
to  a lA%  chord  extension.  The  leading  edge  device  is  12.3%  of  the  extended 
chord  and  is  deflected  60  degrees.  The  main  aft  flap  is  30.7%  of  the  extended 
chord  and  is  deflected  30  degrees.  The  control  flap  was  at  zero  deflection; 
that  is,  the  jet  exited  parallel  to  the  main  aft  flap.  All  the  effects  shown 
are  for  a wing-body  combination  only.  A common  point  of  reference  is  in 
general  given  in  all  the  parametric  effects  studies.  This  common  baseline 
point  is  denoted  by  a circular  symbol.  Some  effects  studies  are  shown  for  off- 
baseline  cases,  for  which  the  corresponding  secondary  variations  are  noted. 


1.  ROLLING  MOMENT  DUE  TO  SIDESLIP 

Figure  35a  presents  the  effect  of  body  size  and  shape  on  the  Incremental 
rolling  moment  derivative,  C*  , duo  to  the  body.  The  effect  of  body  size  is 
shown  for  both  a constant  body  height,  b,  and  body  width,  a.  The  effect  of 
body  shape  for  a constant  body  frontal  area,  Sr,  is  also  given.  The  prediction 
by  Equation  5*2.2. 1-c  of  DATC0m(21)  Is  shown  as  the  dashed  line. 

The  effect  of  body  shape  and  vertical  wing  location  for  bodies  with 
frontal  areas  equal  to  the  baseline  body  is  presented  in  Figure  35b.  As 
expected,  the  magnitude  of  Increases  with  the  high  wing  location.  The  in- 
cremental effect  of  the  body  increases  as  the  wing  becomes  increasingly 
displaced  from  the  center  of  the  body. 

The  effects  of  sweep  and  vertical  wing  location  are  presented  in  Fiqurcs 
35c  and  3$d  for  the  extremes  of  body  shapes  considered.  Note  that  the  Use- 
line  configuration  having  a circular  body  appears  slightly  off  the  parametric 
line  of  Figure  35c.  As  expected,  wing  sweep  angle  is  a very  significant 
parameter. 
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Geometric  dihedral  angle  is  another  significant  parameter  in  the  evalua- 
tion of  wing- body  rolling  moment  due  to  sideslip.  The  effect  of  geometric 
dihedral  and  vertical  wing  location  is  shown  in  Figures  35e  and  3 5f  for 
circular  and  vertical  elliptic  bodies,  respectively.  As  can  be  seen,  a few 
degrees  of  geometric  dihedral  can  offset  the  effect  of  vertical  location  of 
the  wing.  Thus,  for  swept  wing  aircraft  high  wing  configuration  have  negative 
dihedral  or  anhedral  and  low  wing  configurations  have  dihedral  to  maintain 
simi  lar  levels  of  Cj,g. 

The  effect  of  momentum  distribution  is  also  a very  significant  parameter 
in  predicting  the  wing-body  Cj,g.  Previous  difficulty  in  the  analysis  of  this 
parameter  is,  of  course,  the  primary  reason  for  the  development  of  the  present 
method.  In  illustrating  the  effects  of  momentum  distribution  the  fraction  of 
blown  span  and  momentum  injection  angle  were  considered.  In  all  cases  the 
total  wing  momentum  coefficient  was  held  constant  at  Cy  = 1.01.  Figure  35g 
illustrates  these  effects  about  the  baseline  configuration.  Recall  that  the 
baseline  had  a low  level  (section  cy  = 0.263)  of  blowing  outboard  of  the  70% 
span  station.  In  Figure  35g,  the  50%  inboard  span  blowing,  as  well  as  the 
baseline  70%,  has  an  outboard  section  blowing  level  of  cy  = 0.263.  The  jet 
angle,  <5jet,  is  positive  when  the  jet  is  injected  downward  relative  to  the 
flap  reference  plane.  In  this  figure  note  that  the  injection  angle  is  a 
strong  parameter  whereas  the  effect  of  the  percent  span  of  the  main  jet  flap 
is  small,  and  in  general  the  magnitude  of  decreases  as  the  jet  flap  span 
decreases.  The  cases  illustrated  in  Figure  35h  differ  only  in  that  there  is 
no  blowing  on  the  outboard  section  (outboard  potential  flow  assumed).  Note 
that  the  effect  of  jet  flap  span  is  much  greater  and  the  trends  are  opposite 
of  those  shown  in  the  previous  figure.  The  magnitude  of  Cfln  in  this  case 
shows  significant  increases  as  jet  flap  span  is  reduced.  These  two  cases 
presented  in  Figures  35g  and  35h  graphically  illustrate  the  sensitive  nature 
of  the  rolling  moment  derivative  to  small  changes  in  the  distribution  of 
blowing  across  the  span.  Contrasted  to  the  high  wir.q  location  data  of  Figure 
35h,  data  for  a low  wing  configuration  are  presented  in  Figure  35 • • Note  the 
vertical  wing  location  effect  is  nearly  a constant  (AC j u = + 0.0023)  for  all 
combinations  shown.  This  increment  .ilue  is  also  approximately  the  same  as 
shown  in  Figure  35b,  which  indicates  that  body  effects  can  be  superimposed 
and  maintain  a good  degree  of  accuracy. 

The  final  effects  shown  are  for  wing  aspect  ratio  and  taper  ratio  varia- 
tions about  the  baseline  planform.  These  effects  are  presented  in  Figure  35j. 
The  momentum  distribution  for  these  cases  differ  from  the  baseline  In  that 
there  is  no  blowing  outboard  of  the  70%  span  station.  The  baseline  point  is 
denoted  by  the  circular  symbol  for  a point  of  reference.  Note  that  the  taper 
ratio  has  little  effect  on  rolling  moment  derivative  for  this  distribution  of 
momentum.  Aspect  ratio  has  a relatively  small  effect  as  compared  to  other 
parameters  examined. 


2.  SI  DEWASH  DUE  TO  SIDESLIP 


The  parametric  effects  of  these  variables  on  sidewash  derivative  at  a 
typical  field  location  are  presented  in  Figure  36.  The  field  point  selected 
corresponded  to  the  quarter  chord  of  the  vertical  tail  , mean  aerodynamic 
chord  on  the  Ames  I8F  model  described  in  Reference  16.  This  point  was  main- 
tained constant  with  respect  to  the  wing  apex.  Effect  studies  where  wing 
vertical  location  is  a parameter  may  therefore  require  special  interpretation. 
The  sidewash,  as  predicted  by  the  computer  program,  is  based  on  a flat  wake 
trailing  back  streamwise  from  the  one-quarter  chord  of  the  wing.  This  assump- 
tion is  a first  order  approximation  to  the  physical  wake  location  and  in 
general  is  sufficient  for  accurate  downwash  prediction  in  the  vicinity  of  the 
horizontal  tail.  In  the  case  of  sidewash,  however,  the  vertical  location  of 
the  field  point  in  question  relative  to  the  wake  is  a very  strong  parameter. 
Therefore,  some  judgment  in  the  use  of  the  sidewash  prediction  methods  may  he 
requi red. 

The  effects  of  body  shape,  holding  body  frontal  area  equal  to  the  baseline 
value,  and  vertical  location  of  the  wing  are  presented  in  Figure  36a.  Note 
that  the  high  wing  location  shows  a destabilizing  effect  of  approximately  35 
percent.  The  base  case  is  again  denoted  by  the  circular  symbol.  The  effects 
of  wing  sweep  and  vertical  location  are  illustrated  in  Figure  36b.  These 
cases  differ  from  the  baseline  in  that  the  body  shape  is  that  of  a flat 
(b/a  * 0.5)  ellipse.  The  baseline  point  is,  therefore,  off  the  line  of  data 
presented  but  is  given  for  reference.  As  before  the  high  wing  destabilizing 
effect  is  approximately  the  same,  indicating  that  the  effect  of  vertical  wing 
location  could  be  superimposed.  As  expected.  Increased  wing  sweep  has  a de- 
stabi 1 iztng  effect. 

The  effect  of  wing  dihedral  and  vertical  location  is  presented  in  Figure 
36c.  Again,  superposition  of  vertical  location  of  the  wing  appears  to  be 
possible.  Positive  dihedral  increments  would  tend  to  cause  slight  decreases 
in  vertical  tail  stability  contribution  to  restoring  yawing  moment. 

The  effect  of  momentum  distribution  via  jet  Injection  angle  and  percent 
span  of  the  jet  flap  are  shown  in  Figures  36d,  3&c  and  3&f.  Perturbations 
about  the  base  case  are  given  in  Figure  36d,  For  this  type  of  distribution 
the  half  span  Jet  flap  case  shows  a decrease  In  stability  as  Injection  angle 
is  increased.  When  the  small  amount  of  outboard  blowing  present  on  the  base- 
line Is  removed,  the  trend  reverses,  as  seen  in  Figure  36c.  The  superposition 
of  effect  of  vertical  wing  location  is  again  verified  between  the  high  wing 
cases  of  Figure  36e  and  the  low  wing  cases  of  Figure  36f. 

The  effect  of  planform  on  sidewash  derivative  Is  Illustrated  in  Figure 
36g.  These  cases  differ  from  the  baseline  in  that  there  is  no  blowing  out* 
board  of  the  70%  span  station.  The  base  case  Is  denoted  on  the  figure  by  the 
circular  symbol  for  reference  purposes.  The  figure  indicates  that  Increases 
in  aspect  ratio  as  well  as  taper  ratio  would  reduce  the  vertical  tail 
stabilizing  effect. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


A simplified  method  has  been  devised  for  calculation  of  lift,  axial  force, 
rolling  moment,  sldewash,  and  downwash  for  jet  flapped  wing-body  combinations 
In  pitch  and  sideslip.  A number  of  comparisons  have  been  made  with  experi- 
mental data.  From  these  comparisons  the  following  general  conclusions  have 
been  drawn. 

(1)  Lift  coefficient,  as  a function  of  angle  of  attack,  is  predicted 
reasonably  well  over  the  range  where  linearized  potential  theory  may  be 
expected  to  apply.  In  many  cases  this  is  true  well  beyond  the  normal  bounds 
of  applicability.  The  exceptions  are  normal,  such  as  upper  surface  stall, 
undersurface  separation  from  large  leading  edge  flaps,  extremely  high  control 
flap  deflections,  and  low  Cp  Insufficient  to  establish  flow  attachment.  The 
preceding  statements  are  true  over  a wide  range  of  planform  shape,  leading 
and  trailing  edge  mechanical  flap  geometry,  deflection,  and  spanwise  arrange- 
ment, and  level  and  trailing  edge  distribution  of  the  injected  momentum  vector. 

(2)  Conclusion  (1)  generally  applies  to  axial  force  calculation,  although 
it  Is  noted  that  the  planform  variations  in  the  data  comparisons  are  more 
restrictive.  The  correct  assessment  of  axial  force  requires  consideration  of 
the  ultimate  wake  turning  as  well  as  the  spanwise  distribution  of  vortex 
loading.  For  the  cases  cited,  good  results  are  obtained  without  direct  con- 
sideration of  nonplanar  effects  in  the  near  field. 

(3)  Rolling  moments  due  to  sideslip  are  adequately  predicted  for  wing- 
body- vert i cal  tall  combinations  when  wing-body  computer  calculations  are 
corrected  for  the  vertical  tail  contribution.  The  rolling  morr«ent  calculation 
does  not  include  the  effects  of  pylons  and  nacelles,  which  may  have  a 
significant  contribution  on  some  configurations. 

(I»)  The  entire  EMF  concept  appears  to  provide  a good  preliminary  design 
working  tool.  In  view  of  the  general  levels  of  correspondence  for  lift,  axial 
force,  and  rolling  moment,  it  is  believed  that  span  load  distributions  are 
necessarily  correct  to  about  the  same  level. 

« 

Based  on  the  preceding  conclusions  the  following  areas  are  recommended 
for  follow-on  studies. 

(1)  The  input-output  routines  should  be  modified  to  Include  nonsymmetric 
momentum  Injection  and  flap-aileron  deflection.  This  should  provide  Insight 
into  engine  or  duct  failure  and  lateral  trim  against  such  failures.  For  such 
cases  the  yawing  moments  due  to  induced  drag  asymmetry  can  become  important, 
and  can  perhaps  be  reasonably  estimated  by  this  method.  Similarly,  the  im- 
position of  an  antisymmetric  twist  angle  should  give  insight  into  the  roll 
damping  for  jet  flap  type  high  lift  systems. 
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(2)  Pitching  moment,  not  considered  in  the  present  study,  is  also  a 
logical  extension.  The  expected  accuracy  Is  not  clear  because  of  the  simpli- 
fying assumptions;  however,  it  is  expected  that  at  least  first  order  effects 
can  be  calculated. 
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Figure  1.  Essential  Program  Elements 
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Spence's  2-D  Jet  Flap  Model 
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Figure  3.  Coordinates  for  Airfoil  with  Simple  Hinged  Flap 
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EMF  LOAD  » j sin  (o  + 6f  + <5 j ) 


Figure  k.  Schematic  of  2-0  Airfoil  with  Nose  and  Main  Flap 
and  Equivalent  Mechanical  Flap 
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(a)  c'  ts  sensitive  to  effective  a;  location  of  3 A c'  control  point  is 
variable  during  span  load  iteration. 
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Figure  6.  illustration  of  Boundary  Conditions 
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Figure  8.  Pitch  Plane  and  Simple  Sweep  Plane  Resolution  for  Momentum 
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WING  CIRCULATION  LIFT 
EMF  CIRCULATION  LIFT 


o:  = RESIDUAL  LIFT 


= DIRECT  THRUST  LIFT  c^SINe 


SPAN  STATION  n 

(a)  Near  Field  Calculation,  S = aj 


A = WING  CIRCULATION  LIFT 
B = TOTAL  JET  SHEET  CIRCULATION  LIFT 


1 , 

B 
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C = cy  c * 
P avg 


RESIDUAL  LIFT 


B+C  = DIRECT  THRUST  LIFT  cosine 


SPAN  STATION  n 
(b)  Far  Field,  B = Bo, 


Figure  10,  Continued  Conversion  of  Jet  Momentum  to  Vortex  Forces 
Beyond  Near  Field  Lift  Calculation 
(a)  and  (b)  have  equal  lift,  unequal  drag 


Figure  11.  Functional  Flow  Chart 
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Figure  \k.  Comparison  of  Lift  Curve  Slope 
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ANGLE  OF  ATTACK,  a 


Figure  19.  Lift  Comparison,  Large  Scale  IBF,  Full  Span  Jet  Flap  Stowing, 
Landing  Flaps 
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Figure  20.  lift  Comparison,  Large  Scale  I8F,  Fart  Span  Jet  Flap  Blowing, 
Landing  Flaps 
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AR  - 8,  A - 27.5°,  X - 0.3 
«f  - 60°,  6Ail  - 6o°,  6fa  - o° 

FULL  SPAN  JET  FLAP  BLOWING 


MOMENTUM  COEFFICIENT,  Cu 


figure  21.  Lift  Comparisons  at  a • 0,  Large  Scale  ISF 
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Figure  23.  Lift  Comparison,  Large  Scale  IBF,  Part  Span  Jet  Flap 
Stowing,  60°  Main  Flap  Autk.  40°  Control  Flap 
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gurs  25-  General  Arrangement  of  NASA-Ames  EBF  Model 


ANGLE  OF  ATTACK,  DEG.  AXIAL  FORCE  COEFFICIENT 


ANGLE  OF  ATTACK,  OEG.  AXIAL  FORCE  COEFFICIENT 

Figure  27-  Concluded. 
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TABLE  I SUMMARY  OF  SIDESLIP  TEST  RUNS , 
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Figure  32.  Continued 

(c)  Cu  = 1.52,  6f  - 30°,  <$a 1 1 “ 30°,  5fa  « 0° 
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Figure  32.  Continued 

(f)  Cv  = 1.05,  6f  - 60°,  6a i i - 30°,  6fu  " 20° 


a = 4° 
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WING/BODY  6 ESTIMATED 

VERTICAL  TAIL 
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Figure  32.  Concluded 

(h)  Cu  - 1.03,  6f  - 60°,  Sail  B 30°,  6fu  * 
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Figure  33.  Rolling  Moment  Comparison,  Large  Scale  USB 
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Figure  35.  Continued 

(h)  Effect  of  Manehtum  Distribution;  Baseline  with  NO 
Outboard  Blowing 
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Figure  35*  Concluded 

(1)  Effect  of  Planform  Variations 
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Figure  36.  Effect  Studies  - Sldewash  Derivatives 
(a)  Effect  of  Body  Shape  and  Location 
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Figure  36.  Continued 


(e)  Effect  of  Nooentua  Distribution;  Baseline  with  No 
Outboard  Blowing 


S I DEWASH 
DERIVATIVE 
do/ dS 


Figure  36.  Concluded 

(g)  Effect  of  Planfonn  VarUtJos;  Baseline  with  Mo  Outboard 


APPENDIX  A 


CALCULATION  OF  BODY- INDUCED  ANGLE  OF 
ATTACK  ON  THE  WING 


The  geometries  of  the  elliptic  cross  section  and  its  circular  transforma- 
tion are  shown  in  Figure  A-l.  The  sideslip  crossflow  velocity  VsinB  a o the  (v,w) 
and  (v',w')  velocities  are  also  noted.  The  case  shown  is  for  the  long  dimen- 
sion vertical  (b/a  > 1),  for  which 


transforms  the  ellipse  to  a circle.  The  complex  potential  in  Mesllp  is 


F(c)  = -V  sin8[c  + ~ j 

The  complex  velocity  in  the  real  plane 
direction) 

dF  . dF  1 

3z  ’ v * ,w  * 3?  dF" 
3s 

From  Egs  A. 2 and  A. 3-  we  obtains 


is  (note  that  VsinS  is  in  negative  y 


dF 

dc 


SSLmJLm  ISl*.  

s V®  (C2  - n2  ♦ a2r  * iJV 

where 
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and  where 


= F 


C2  + n 


a2ri 

■ , z = n + , " 

2 £2  + n2 


A. 6 


The  induced  angle  Aa/f5  in  the  real  plane  is  written  in  transformed  coordinates 
(n,S)  since  a quart ic  solution  is  required  to  express  it  in  real  coordinates. 
However,  the  transformed  coordinates  (n»5)  are  quickly  obtained  by  writing 
Eq  A. 6 as 


(*-)  + 
al 


(f)  (f) 

al  al  i-1 


L.V2  + (1L-)2 
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’1  i-1 
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ai  al  1-1 


<4r> 2 + (rr> 


al 


i-i 


31  i-1 


A. 7 


and  iterating  simultaneously  in  each  step, starting  with 


The  iterated  results,  with  Eq  A. k , give  the  induced  angle  Aot/g  as  a function 
of  (y,z)  coordinates. 

In  pitch,  the  crossflow  velocity  is  rotated  90  degrees  from  the  sideslip 
case, for  which  the  complex  potential  is 


F(c)  ■ 


i V sinabody 


9; 


t 


With  Eq  A. 3;  , but  with  dF/dZ  v - ?w,  tnis  gives 


+ (E2  -ti2)  (R2+  a2)  + R2a2 


(E2  ~ n2  + a2)  + 4 52n2 


where  Eqs  A. 5 » A. 6 , and  A. 7 apply. 

For  bi/ai  < 1,  the  appropriate  transformation  is 


Z » c + - 
C 


However,  since  a2  is  replaced  by  -a2, it  also  follows  that 


(r-)2  = 


<V  - i 

'ar 


so  that  algebraic  signs  compencate  when  the  preceding  results  are  used  intact. 


APPENDIX  B 


CALCULATION  OF  AXIAL  FORCE  COEFFICIENT 


The  far  field  wake  Is  assumed  to  consist  of  a rolled  up  vortex  system 
which  Is  coaxial  with  the  redirected  jet,  deflected  downward  by  an  angle  £«,. 

The  vortex  strength  on  each  side  is  denoted  by  r0,  where  T0  is  the  centerline 
circulation  around  both  wing  and  sheet,  and  is  therefore  the  wake  vorticity 
accumulated  by  shedding.  The  total  jet  momentum  in  the  far  field  Is  assumed 
to  be  nCy,  where  ri  Is  the  usual  "turning  efficiency." 

The  geometry  is  illustrated  in  Figure  B-1.  The  jet  and  vortex  are 
approximated  by  a line  vortex  and  jet  as  shown,  and  the  two  sides  are  separated 
by  a distance  (it A) bAi,  where  Ai  is  to  be  determined.  The  vortices  are  in 
equilibrium  at  an  angle  g®  when  the  freestream  crossflow  velocity  balances  the 
mutually  induced  velocity,  giving 


r 

sing  * K 2L 

“ 2bVAi 


B.l 


where  K = b/v1  for  the  rolled  wake  shown,  and  where  K * 2 for  a flat  wake. 

The  vortex  o - circulation  lift  is  the  transverse  momentum  per  unit  length 
of  vortex,  and  is  given  bv 


'•r  J ?Vro  f Alb 

from  which 


ro 

JbVA, 


B.3 


Integration  of  lifting  vortex  loads  directly  gives 


Lr  «•  pV"  o 


B.4 


where  f is  the  average  T across  the  span.  From  Eqs  B.2  and  B.A  it  Is  seen 
that  Ai  * kh  f/r0.  It  is  noted  th..:  ■ I for  an  elliptic  vortex  distribu- 

tion. 
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The  vortex  "drag"  force,  Inclined  by  3 to  the  freestream,  can  be  shown 

to  be 


Cpyortex  = ,po_  2 
tt/JR  B 2bV 


B.5 


where  B Is  the  ratio  of  vortex  drag  to  that  for  an  elliptic  vortex  distribution 
at  the  same  To/2bV.  Substitution  from  F.q.  c.3  allows  the  preceding  to  be 
written 


where  ej,  is  the  span  efficiency  of  the  total  vortex  system. 

With  this  In  mind,  the  total  vorticlty  in  the  wake  is  given  by 

r°Total  _ Clr 
2bVA  " tt/R  A*2 

C^Totai  _ nCu  . 
n/R  Aj  *-'  *AR  Aj  - S ° 00 


- ClTotaI  . nCu  ^ rQTotal 
tt/RAj2  2bVAl 

from  which 
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B.7 
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The  pitch  plane  projection  of  the  injected  jet  is  denoted  by  riCy  In 
keeping  with  the  usual  efficiency-turning  polar  notation.  The  Trefftz  plane 
efflux  of  pitch  plane  momentum  is  assumed  to  be  nC  and  the  lateral  compon- 
ents are  ignored.  This  can  be  done  since  the  effeHtive  "span"  of  the  lateral 
components  Is,  in  the  present  context,  the  vertical  dimension  or  thickness  of 
the  jet.  This  corresponds  to  an  extremely  high  value  of  jj  /[q  (thickness)2) 
(analogous  to  very  high  C /A1 ) , so  that  substantially  no  bend-back  and  assoc- 
iated thrust  recovery  and^vortex  drag  are  expected.  The  axial  force 
coefficient  is  written  (see  Figure  10) 


nAl  A^2 


jr/R  A, 2 


uvortex 

w/RAj2 


which,  for  small  6a,  becomes 


Cx  nC  fnC  g 2 Cq 
* a l1  - u + ^vorl 

w/RAj2  ir/RA2  ir/RAj"  2 ir/RA2 


('  -%-) 


From  previous  relationships  (see  Eqs  B.1,  B.5,  and  B.6),  it  is  seen  that 
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Substitution  of  these  relationships,  with  Eq.  8.7  , into  Eq.  B.8  gives, 
after  some  manipulation,  . 
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Terms  of  order  $ and  C /ufl*  are  usually  small  compared  to  unity.  By  neg- 
lecting their  squares  land  products  and  accounting  for  profile  and  ram  drag 
the  axial  force  expression  B.9  becomes 


1 


nC„  - Cr>  - cD 
» dram  p 


uAt  St 


P 


2KnC  , 

+^(1 


K er, . 


B.  1 U 


In  two  dimensional  flow  the  axial  force  expression  B.10 

corresponds  to  a 100  percent  recovery  of  the  effective  gross  thrust  nCu, 
which  is  reduced  by  ram  and  profile  drags.  Accordingly,  the  last  term  is  a 
total  drag  due  to-  the  three  dimensionality  of  the  flow,  and  Is  therefore 
properly  considered  as  Induced  drag.  The  induced  drag  consists  of  the 
induced  drag  in  the  normal  sense  (extended  to  include  the  effects  of  jet 
sheet  vortex  system)  plus  a loss  of  thrust  recovery  resulting  from  failure 
of  the  jet  efflux  momentum  to  ultimately  align  Itself  with  the  frcestream 
direction. 

It  will  be  observed  that  the  induced  drag  at  Cu  = 0 corresponds  to  the 
usual  expression 


CD-  “ CL2/**er  B.11 


The  span  efficiency,  er*  is  of  course  associated  with  the  spanwlse  distri- 
bution of  circulation  lift  over  the  combined  wing  and  vortex  sheet.  How- 
ever, an  additional  distribution  factor  Aj  is  associated  with  tie  C term. 
The  presence  of  two  parameters  ej,  and  Aj  in  the  Induced  drag  expression 
(in  contrast  to  the  usual  single  parameter  e^)  is  explained  as  follows. 

The  classic  expression  for  induced  drag,  Equation  B.11,  is  based  on  a 
small  wake  Inclination  B.  When  rigorously  developed  from  Trefftz  plane 
considerations,  this  result  is  in  fact  independent  of  the  wake  angle,  B, 
parameter.  This  is  true  for  the  Induced  drag  of  vortex  systems  associated 
with  either  a wing  or  a wing  plus  jet  sheet.  The  additional  induced  drag, 
which  is  the  three  dimensional  loss  of  thrust  recovery,  ts 


It  is  directly  dependent  upon  B,  where  B is  partially  defined  by  the  term 
A)  (Eq.  3.1).  It  ts  noted  that  Aj  is  the  spanwlse  distance  between  "cen- 
troids" of  the  left  and  right  hand  trailing  vortex  systems,  normalized  on 
the  corresponding  distance  for  an  elliptic  bound  vortex  distribution.  In 
the  literature  It  is  generally  referred  to  as  the  "centroid  of  impulse". 


V 
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Equation  B.10  was  developed  on  the  assumption  of  a fully  rolled  wake 
for  which  K = 4/ir2.  However,  if  the  wake  Is  takei.  to  be  flat,  K * 2,  in 
which  case  Eq.  B.10  is  still  valid.  For  a flat  wake  and  elliptic  loading 
(fhi  ® 1,  er  » 1),  Eq.  B.10  reduces  to  the  Maskel  1 -Spence  equation.  As  long 
as 


/ 


the  axial  force  is  relatively  insensitive  to  K.  For  large  C , with  the 
blowing  concentrated  over  a portion  of  the  span,  this  is  not^true,  so  a 
realistic  value  of  K must  be  used.  It  is  believed  that  K = k/n2,  for  the 
rolled  up  wake,  is  most  representative  and  is  the  basis  for  all  axial  force 
calculations  presented  in  the  main  text. 

Utilization  of  Eq.  B.10  requires  that  Aj  and  ef  be  calculated.  It  is 
emphasized  that  these  are  distributional  quantities  associated  with  the 
total  vortex  system.  Some  approximation  is  therefore  necessary. 

The  span  load  distributions  calculated  by  the  method  presented  in  the 
main  body  of  the  report  consist  of  a vortex  loading  on  the  wing  and  EMF 
plus  a residual  lift  term,  corresponding  to  K * 1 , for  which  the  conversion 
of  jet  momentum  to  sheet  vortex  loading  is  incomplete.  The  appropriate 
vortex  distribution  for  calculation  of  A*  and  er  is  therefore  somewhere 
between  the  vortex  lifts  and  total  lift  distributions  as  presently  calcu- 
lated, and  illustrated  in  Figures  (10a)  and  (10b)  in  the  main  text. 

The  differences  in  Aj  and  e evaluated  for  these  extremes  are  typically 
small;  both  quantities  are  therefore  based  on  the  vortex  lift  distribution 
with  K =>  1. 
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Figure  C-1.  Far  Field 
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for  Axial  Force  Calculation 


APPENDIX  C 


COMPUTER  PROGRAM  INPUT,  OUTPUT,  AND  EXAMPLE  PROBLEM 


1.  INPUT  FORMAT  AND  SPECIAL  INSTRUCTIONS 
a.  Basic  Input 

The  total  number  of  Input  cards  varies  according  to  both  the  number 
of  points  used  tc  describe  the  platform  and  the  number  of  field  points  at 
which  flow  angularity  is  to  be  cahulated.  Card  1 is  a title  card  containing 
70  columns  of  free-form  alphameric  data  which  is  printed  out  as  the  output 
head i ng . 

The  remaining  card  input  is  numerical  data.  The  majority  of  these  inputs 
are  in  a floating  point  format  of  seven,  10-column  fields  per  card  (7F10.0) 
and  are  indicated  by  "Field  1",  "Field  2",  etc.,  in  the  following  instructions. 
The  remaining  input  is  in  integer  form  (212)  and  Is  indicated  in  the  following 
instructions  by  "Columns  1 and  2",  etc.  All  integer  input  is  "right  justified", 
and  must  be  located  as  far  as  possible  to  the  right-hand  side  of  the  allocated 
columns. 

The  input  for  each  card  follows.  In  certain  places  reference  is  made  to 
Note  1,  Note2,  etc.,  which  follow  the  card  description. 


Card  1:  Title  (70  column  alphameric  input). 


Card  2: 


Columns  1 6 2;  NY,  number  of  span  stations  at  which  section 

properties  will  be  Input  (right  justified).  This  is  for 
one  wing  panel  only,  and  will  include  the  theoretical 
wing  root  and  tip.  2 5 NY  < 20. 

The  following  block  of  cards  contains  section  characteristics  for  the  NY 
number  of  points  specified  in  Card  2.  Two  cards  are  required  for  each  point. 
Cards  3 and  k are  for  the  wing  centerline.  Subsequent  cards  are  input  in 
order  of  spanwise  location,  ending  with  the  wing  tip.  The  locations  of  the 
Input  points,  while  arbitrary,  should  be  chosen  to  include  those  span  stations 
at  which  section  characteristics  change  rapidly  or  are  discontinuous.  At  such 
locations,  properties  should  be  input  both  just  inboard  and  just  outboard  of 
discontinuities  In  geometry  and/or  at  blowing  levels.  Section  properties  at 
non-input  locations  are  obtained  by  linear  interpolation.  All  inputs  used  are 
for  mechanical  f laps  extended  but  rotated  back  into  the  wing  plane  and  are  In 
£ stream* I sc  reTe fence  sys tern. 
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Card  3:  Input  the  following  for  wing  centerline: 

Field  1 dimensional  value  of  y coordinate 

2 dimensional  cnord,  all  flaps  extended 

3 section  incidence  with  respect  to  the  fuselage  reference 
(in  degrees) 

4 section  zero  lift  angle  of  attack,  in  degrees  (see  Note  1) 

5 main  flap  chord:  wing  chord  ratio  (ail  flaps  extended) 

6 section  momentum  coefficient  (see  Note  2) 

7 main  flap  deflection,  degrees,  positive  trailing  edge  down 


Card  4:  Field  1 6j , jet  momentum  injection  angle,  degrees,  positive  down- 

ward from  flap 

2 e,  jet  skew  angle  (see  Note  2) 

3 C^o,  section  lift  coefficient,  unflaped,  at  zero 
angle  of  attack 


4 leading  edge  flap  chord:  wing  chord  ratio 

5 leading  edge  flap  deflection,  degrees,  positive  nose  down 


Card  5 , 
Card  6: * 


Same  as  for  Cards  3 and  4 but  for  second  point  of  definition. 


Card  8'^  R*P®at  preceding  input  for  third  point  of  definition. 

Card  2 (NY)+1  . „ ..  , , . 

Card  2(MY)+2  ReP«a*  preceding  input  for  wing  tip. 

Card  2(NY)+3: 

Field  I body  width  or  diameter,  dimensional 

2 wing  height  above  body  center,  dimensional 

3 body  height  for  elliptic  bodies,  dimensional,  input  blank 
or  zero  for  circular  body 
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t * 


Card  2(NY)+A 


Field  1 drag  coefficient;  total  of  configuration  profile  drag  plus 
ram  drag.  (Used  in  final  axial  force  calculation)  Based 
on  reference  area. 

2 turning  efficiency,  n>  rejection  of  resultant  efflux 
momentum  into  pitch  plane  and  normalized  on  input  total. 

3 increment  of  semispan  at  which  loading  will  be  output, 
i.e.,  0.05,  0.10,  etc. 

k tOKPCT,  permissible  percent  RMS  error  in  calculated 
effective  angle  of  attack  for  all  control  points  (see 
Note  l),  suggested  value  is  2.0  percent. 


Card  2 (NY) +5: 

Columns  1 & 2;  NC,  number  of  control  points  per  semi  span  used  in 
calculation  procedure.  (Right  justified).  This  number 
will  include  a centerline  control  point  and  a tip  control 
point.  4 < NC  < 16.  Total  used  in  sideslip  calculation 
Is  2 (NC) -1 . (see  Note  3) 

Columns  3 S **;  NPTS,  Number  of  field  points  where  sidewash  and 
downwash  is  desired  (Right  justified).  0 < NPTS  < 20, 
blank  or  zero  if  none  required. 


Card  2(NY)+6:  Lccation  of  first  sidewash/dr.nwash  field  points,  dimensional, 

optional  if  blank  or  zero  input  in  Columns  3_i»  of  previous 
card.  All  locations  relative  to  wing  apex. 


Field  1 y 

2 Y 

3 z 


(positive  towards  aft  fuselage) 
(positive  out  R.H.  wing) 
(positive  up) 


Card  2(NY)+7: 


x.y.z  location  of  second  field  point. 


Card  2{NY)+5+NPTS  Location  of  last  field  point. 
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Card  2(NY)+6+NPTS 

Field  1 wing  span,  dimensional 

2 aspect  ratio  (based  on  span  and  desired  reference  area) 

3 wing  sweep  (degrees) 

4 dihedral  (degrees) 

5 flap  hinge  line  sweep  (degrees) 


Card  2(NY)+7+NPTS 

Field  1 ALPHAS,  angle  of  attack  of  the  fuselage  reference  line, 
degrees 

2 sideslip  angle,  degrees,  positive  far  incident  wind  from 
pilot's  right 

3 FACTOR,  convergence  factor.  Normally  leave  blank,  otherwise 
see  Note  1 . 


b.  Special  Notes. 

(1)  Note  1:  Input  Parameters  to  Aid  in  Program  Convergence.  For 

most  normal  caseruns,  the  three  input  parameters  which  deal  with  convergence 
are  input  as  zero  or  blank. 

Lards  3,  4.  ...2(NY)+2,  Field  4 

Section  angle  of  attack  for  zero  lift  (flapped  and  powered).  This 
Input  parameter  Is  used  only  for  setting  up  an  initial  point  for 
the  convergence  procedure.  Normally,  zero  values  for  this  parameter 
are  sufficiently  close  for  convergence.  For  cases  where  convergence 
is  not  obtained  (output  ESTAR»E0K) , a more  accurate  estimate  of 
this  parameter  may  be  required.  Variations  of  this  parameter  is  the 
most  effective  means  of  solving  convergence  problems. 


Card  2(NY)+4,  Field  4:  permissible  error,  EOKPCT 

If  left  blank  or  Input  as  zero,  the  program  internally  sets  the 
value  at  1|  of  the  input  angle  of  attack  of  the  fuselage  reference 
line,  EOK  » (J .}( .01 ) (ALPHA?) . For  many  cases  a larger  permissible 
error  may  yield  equally  valid  results  thus  requiring  less  iterations 
and  computer  time. 
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Card  2(NY)+7+NPTS,  Field  FACTOR,  convergence  factor. 


Normally  left  blank  or  zero,  whereby  the  computer  Internally  sets 
FACTOR  to  0.4.  If  the  first  attempt  at  convergence  falls,  the 
computer  changes  FACTOR  to  0.8  and  starts  over.  If  convergence  Is 
not  obtained  on  the  second  time  (finat  output  of  ESTAR>E0K) , the 
case  may  be  rerun  with  another  Input  value  for  FACTOR.  Recommended 
range  is  greater  than  0.2  but  less  than  2.0.  Note:  Do  not  set 

FACTOR  equal  to  1.0. 

The  n1^  term  of  the  nc^  guess  is  multiplied  by  FACTOR.  Do  not  set  FACTOR 
equal  to  1.0  since  a singular  matrix  results. 

(2)  Note  2:  Momentum  Coefficient.  It  is  important  that  the 

trailing  edge  momentum  vector  be  properly  defined,  since  (a)  its  projection 
into  the  simple  sweep  plane  is  required  for  iterative  determination  of  EMF 
size  and  total  circulation  and  (b)  its  projection  into  the  pitch  plane  is 
required  for  calculation  of  residual  lift  and  for  calculation  of  axial  force 
coefficient.  The  momentum  coefficient  given  in  Card  3»  Field  6,  is  the 
magn'tude  of  the  local  momentum  efflux  vector  Cy  per  unit  of  physical  span  (y 
direction)  which  leaves  the  trailing  edge,  normalized  on  freestream  dynamic 
pressure  and  freestream  chord.  This  is  the  total  per  unit  of  span  regardless 
of  the  direction  of  the  efflux.  The  *et  skew  angle  e of  srd  4,  Field  2 is 
the  angle  between  the  local  efflux  vc  or  and  the  simple  sweep  plane. 

The  Input  angles  e,  Sf,  6 \ , and  Ac/2*  are  sufficient  for  the  required 
resolution.  For  IBF  systems  the  eft  lux  from  slot  blowing  is  substantially  in 
the  normal  plane;  for  such  cases  one  should  input  e as  zero.  In  this  case  the 
pitch  plane  projection  or  "turning  efficiency"  n Is  a fallout  quantity 
(scrubbing  losses  are  assumed  to  be  accounted  for  in  the  input  Cy) . 

For  cases  other  than  the  18F  the  skew  angle  e is  a spanwise  variable  which 
is  generally  not  well  known.  The  required  assumptions  concerning  efflux 
spreading  are  left  to  the  user's  experience,  his  available  flow  visualization 
and  static  force  data,  and  his  knowledge  of  wind-on  interactions  with  static 
spreading  and  turning.  In  general,  the  input  t and  Cy  should,  when  Integrated 
over  the  sp^n,  be  consistent  with  the  iny«;  "turning  efficiency"  n. 


(3)  Note  3:  Spanwise  Control  Points;  Number  and  location.  The 

spanwise  location  of  the  control  points  are  specified  for  each  number  of 
control  points  chosen.  These  are  tabulated  for  quick  reference  in  Table  C-J 
for  four  through  sixteen  control  points. 

Host  powered  high  lift  systems  are  characterized  by  sharp  breaks  in  flap 
chord,  flap  deflection  and  momentum  distribution.  Control  point  locations 
Should  be  chosen,  wherever  possible,  to  give  a more  or  less  balanced  recogni- 
tion to  the  various  segments  between  these  breaks.  Best  results  have  been 
obtained  when  major  segments  of  the  span  have  approximately  the  same  number 
of  control  points.  It  is  also  desirable  to  avoid  having  control  points 
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located  too  near  sharp  breaks.  In  certain  cases  it  may  be  desirable  to  run 
more  than  one  number  of  control  points  and  examine  the  resulting  span  load 
distribution. 

2,  OUTPUT  FORMAT 


Line  1 


TITLE  as 

input 

Lines  2 & 3 

ALPHA 

- 

Reference  angle  of  attack,  degrees. 

BETA 

- 

Sideslip  angle,  degrees. 

Sweep 

- 

Sweep  angle  of  the  wing  one-quarter  chord  - degrees 

AR 

- 

Wing  aspect  ratio. 

Span 

- 

Wing  span,  dimensional  as  input. 

Lines  /(  & 5 

CL 

- 

Wing  lift  coefficient  in  the  presence  of  the  body. 

CL(CIRC) 

- 

Wing  circulation  lift  coefficient. 

CDP 

- 

Wing/body  profile  drag  coefficient,  as  input. 

CX 

- 

Wing  axial  force  coefficient. 

ETA 

Turning  efficiency,  q,  projection  of  resultant  effli 
momentum  into  pitch  plane  and  normalized  on  input 
total . 

CMU 

- 

Blowing  coefficient,  Cy. 

CR(ROLL) 

- 

Wing  rolling  moment  coefficient. 

Table  1 

ETA 

- 

Wing  station  as  fraction  of  semispan. 

GAMA /V 

- 

Wing  section  circulation  non-dimensional ! zed  on 

freestream  velocity 
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CL*C/CAVG 
ALPHA  E 


Wing  load  distribution  parameter  (section  lift  per 
unit  span  non-dimensional i zed  on  average  wing  chord) 

Section  effective  angle  of  attack,  degrees. 


Table  2 

X POINT,  - x,  y,  z coordinates  of  downwash,  sidewash  points, 
Y POINT,  as  input. 

Z POINT 


EPSILON  - Downwash  angle,  degrees. 

SIGMA  - Sidewash  angle,  degrees. 


3.  SAMPLE  CASE  INPUT  S OUTPUT 

An  initial  set  of  control  cards  are  required  by  the  CDC  6600  computer  at 
Wright-Patterson  Air  Force  Base,  Dayton,  Ohio.  The  card  deck  as  supplied  by 
the  contractor  will  compile  the  source  program  and  save  the  object  deck  under 
the  catalog  name  "HOLMES, CY=V.  Once  this  is  completed  one  must  estimate  the 
computer  run  time  in  seconds  by  the  following  approximate  formula: 


TIME  "*  7 + 10*  Number  of  cases, 

a case  being  each  a, 6 combination.  For  the  sample  case  (one  a,B  combination), 
TIME  s 17sec. 


The  required  control  cards  and  input  data  are  shown  In  Figure  C-l.  The 
"SETCORE."  card  must  be  used  to  initialize  the  core  storage  to  zero.  The 
resulting  output  is  shown  in  Figure  C-2. 


NUMBER  OF  CONTROL  POINTS  PER  SEMI  SPAN 


4 

5 

6 

7 

8 

9 

10 

0 

0 

0 

0 

0 

0 

0 

0.3827 

0.3090 

0.2588 

0.2225 

0.1951 

0.1736 

0.1564 

0.7071 

0.5878 

0.5000 

0.4339 

0.3827 

0.3420 

0.3090 

0.9239 

0.8090 

0.7071 

0.6235 

0.5556 

0.5000 

0.4540 

0.9511 

0.8660 

0.7818 

0.7071 

0.6428 

0.5878 

0.9659 

0.9010 

0.8315 

0.7660 

0.7071 

0.9749 

0.9239 

0.8660 

0.8090 

0.9808 

0.9397 

0.8910 

0.9848  0.9511 

0.9877 


TABLE  C-1  LOCATIONS  OF  CONTROL  POINTS  IN 
FRACTIONS  OF  SEMI  SPAN 


(a)  Four  through  Ten  Control  Points  per 
Semi  span 
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NUMBER  OF  CONTROL  POINTS  PER  SEMI  SPAN 


11 

12 

13 

0 

0 

0 

0.1423 

0.1305 

0.1205 

0.2817 

0.2588 

0.2393 

0.4154 

0.3827 

0.3546 

0.5406 

0.5000 

0.4647 

0.6549 

0.6088 

0.5681 

0.7557 

0.7071 

0.6631 

0.8413 

0.7934 

0.7485 

0.9096 

0.8660 

0.8230 

0.9595 

0.3239 

0.8855 

0.9898 

0.9659 

0.9350 

0.9914 

0.9709 

0.9927 

14 

15 

16 

0 

0 

0 

0.1120 

0.1045 

0.0980 

0.2225 

0.2079 

0.1951 

0.3303 

0.3030 

0.2903 

0.4339 

0.4067 

0.3827 

0.5320 

0.5000 

0.4714 

0.6235 

0.5878 

0.5556 

0.7071 

0.6691 

0.6344 

0.7818 

0.7431 

0.7071 

0.8467 

0.8090 

0.7730 

0.9010 

0.8660 

0.8315 

0.9439 

0.9135 

0.8819 

0.9749 

0.9511 

0.9239 

0.9937 

0.9781 

0.9569 

0.9945 

0.9808 

0.9952 

TABLE  C~’  (Concluded) 

(b)  Eleven  through  Sixteen  Control  Points  per  Semispan 
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UNIVAC<T17*CM60000)  ACCOUNT * ENGINEER  t PHONE 

ATTACH ( LGO* HOLMES* CY=4 ) 

SETCORE. 

LGO. 

7/8/9  (OVERPUNCH) 

FOL/ARL  CHECK  CASE  * AMES  U/H  IBF  BASELINE*  DEL  F=30. 
4 


0.0 

9.405 

0 

0.  i 

.307 

1.0294 

28.5 

0.0  ' 

0.  ) 

.338 

.1228 

57.0 

15.0 

4.796 

0.0 

0.0 

.307 

1.0294 

28.5 

0.00 

0.0 

.338 

.1228 

57. 

15.01 

4.796 

0.0 

0.0 

.307 

0.26  8 

28.5 

0.0 

0.0 

.3  8 

.1228 

57.0 

21.4475 

2.82 

0.0 

0.0 

.307 

0.2628 

28.5 

0.00 

0.0 

.338 

.12  >8 

57. 

5.833 

2.91 

5.833 

.25 

.96 

0.2 

2.0 

09  4 
?*.791 

3.769 

9.042 

24.06 

0.0 

9.042 

24.06 

0. 

4.296 

24.06 

0.0 

0. 

42.895 

8.0 

27.5 

0.0 

23.0 

4.0 

99. 

6Z7/8/9/ 

1.0 

(OVERPUNCH) 

Figure  C-l . 

Sample  Problem  Input 

ALPHi=  )i 

:STAv=  . 2 ? 3-  ►'•Of  +•']? 

£ S T A 9 - . o .■>  o - . 3 ..  4 u U 
£ S T A ■•'  = . - .7.  >r  i.-  :? 

I T T t.  •*. - 7 -,0<  = . f 


A •, » ■ l ! * 


' t c •<  = . o n *♦  -i  d 


F n L / A L fh  c<  HA'-' 


J/  = l 1 4F  asal!.c  , Or  L F = i .1 . 


L 

c l ■'  ’ (. 

4 _ ■ 4 

; r 0 •?  fl  K 

4 

...  C 1.0.0 

27.  £ o 

H.oo  -2  . c 9 

l :l  ( 

o i F.c  o ; )'• 

:T  t 0 -H.l 

- Z 2 c.'  l < . 

2-94.;  . 7 ~ 1 0 .» 

• 2£&i  C* 

. 9tu  1.00970 

. T 3 

t )•■•:•./  1 

' * / ’1  V«l  , 

> w .4  » •.  4 X x 

AL'J-A  _ 

-1 . 00 0 0 u 

C • • j J 0 0 .. 

1 • 0 o C u *i 

— 2?»lc>77  i* 

• . {?  0 1 '!  (1 

L .27-2“ 

2.33C?  a 

-.1739? 

-.LOO  0 0 

r . '--I--.? 

1.4DV  1 1 

-3. 0^30? 

- .40100 

1 1'  . 2 ft 4 4 7 

4 . u?  9 1 1 

-4. 19622 

-.20  K'O 

1C . - 72?  1 

— . 5 1 'J  J 

- l . - v 7 r,  0 

0 . t 0 O 0 o 

11.72-M7 

*♦.41  2 ? h 

-c . 2 1221 

. 20  0 0 !» 

10.4  0 Of* 

4 . 4 1 JO  3 

-t,.  7 1192 

.MU')  0 0 

1 C . ~ 7 ' 

4,10217 

-4,39017 

. c i 0 0 0 

9.1)11!' 

5 • 52  ftc  5 

-7.22723 

. S 0 0 0 u 

£.40344 

2.3SS71 

3.99322 

i . io  o Jo  r. 

il  1 ) . J i 

J . J 0 0 0 i) 

-22. 9l)9  > 

X POINT 

Y POI  IT  ? °0  IN  T 

r.  PS  I LON 

SIGUA 

CM -OLD 
- . 0 1625 


2l.7  4 1 ,)  0 

2-.  Or>  )M. 
24  , CoUClG 
?4 » 0-j  ou  n 


3 . 7 C 9 ,)  0 
0 . G L Li  L. '! 
0.  Ci GO nu 
0.00001 


0.0-21) 
C. 0-2 j J 
-.2  9ov.l 
0.00)10 


Hi.  coil  Cl 
in.  eq‘J)2 
1 i . 2 2 7 r.  2 
13. 47049 


2.  1G?1  7 
• G .3  9 •)  C 
— . 3^L> 4 c 
. 9767? 


Figure  C_2.  Sample  Problem  Output 


APPENDIX  D 


COMPUTER  PROGRAM  LISTING 


The  original  computer  program  was  written  In  Foreran  V fo.  the  Uni vac  1106 
DEMAND  System.  The  program  listing  which  follows  hrs  been  adapted  to  the  CDC 
6600  computer  at  Wrlght-Patterson  Air  Force  Base. 
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Preceiii*  pift  link 


PROGRAM  UNIVAC ( INPUT » OUTPUT ► TAPEb= INPUT* I' APE6=0UTPUT ) 
COMMON  NS  PAN*  B t ALPHAA  (3  i)»A(3  >3  )*C  (33). 

1 G(3.)*X(33.  5) .ERROR (3  * 3)*E(3  )*ALPHAE(3  ) 

2 * ALFA  II <3  )) > YTAB ( 20 ) »ALPZL(20) 

3 > NY* CLAM 

b * CMUT AB ! 20 ) *CF0CIN<20) * DELPTB (20 ) »DELJTB(2Q) .KWAKE 

6 *EPSTAB(20) *SWE  P*  LAMHL*  BET  A * SGAMA  * CLOTAti  ( 20  ) 

7 * CNOC IN(20 ) * DELNTB (20 ) 

REAL  I vi  TAB  (20  ) > Iw  * LB  ( 3 i * 33  ) .LAMHL 
INACf-TR  BM*  BMP1  * BMP2 

Si  ME  NS  I ON  CTAO(20)*  ALFAE1  (3  ) * ETA  1 1 ( 3.5 ) * 

1 ETA  (40  > * PHI  I (3  ) *ALPHAI(3.:>)  *GAMA0V(3  )* 

2 CLCO.A(iJ) ♦ CO ICOC ( 3 >) . 

3 ALFA  I ( 3 ) *HSTAR(33.3  ) 

4 * PHI  12(3  ) * PHI ( 40 ) * FBAR ( 3 *;).GBAR(3  .3  ).SUM48(3  ) 

5 * CROLl ( 3 ) * T 1 7LE ( 10 ) 

6 *XPB(20) . YPB(20) *ZPB(20) 

READ  (5*  l:  2)  (TITLE(KK  )*K  ••  =1*8  ) 

READ  (b*4)  NY 
RA0=57.296 
4 FORMAT  (312) 

PRINT=0.0 

10  DO  lb  I = 1 * NY  * 2 

READ  (5. 26)  YTAB  ( I ) *CTAB(  I)  . IviTAB  ( I ) *ALPZL(  I) 

+ »CFOCIN(  I ) * CMUTAB  ( I ) * DELFTI3  ( I ) * DEL JTU  ( I ) 

+ * EPST  AB ( I) .CLOTAB(I) .CNO^INt I) * DELNTB ( I) 

READ  (b  *26  ) YTAB  ( 1 + 1 i *CTAB(  1 + 1 ) * IvJTABt  1 + 1 ) *ALP2L(  1 + 1 ) 

+ *CFOCIN(  1 + 1 ) .CMUTAB  ( 1 + 1 ) . DELFTB  ( 1 + 1 ) . DEL  JTB  ( IH) 

+ * EPST  AB  ( 1 + 1 ) *CLOTAB(  1 + 1 ) .CNOdNf  1 + 1 ) * DELNTB  ( I + 1 ) 
lb  CONTINUE 
20  FORMAT  (bFlO.l) 

2b  FORMAT  (6F1U.ii) 

2o  FORMAT (7F10.  O 
KWAKE=1 

READ  (b*2e>)  DIAB.  ZWING.BMAJOR 
READ (b. 26)  COOT  * ETACMU*  Db'LY  » EOKPCT 
IF ( EOKPCT . LE . 0 i)  EOKPCT=l .0 
READ  ( b *4 } NN» NPTS 
IF(NPTS.LE.U)  GO  TO  3 4 

DO  3 KP=1»NPTS 

READ (b *20 ) XRB(KR) *YPU(KP) *ZPH(KP) 

3 CONTINUE 

3 34  HEAD  t b * 2 U ) B.  AH.  S/iE  P*GAMA, LAMHL 
T APEHrC  TAB  ( NY ) /CTAfi  ( 1 ) 

SG  AM  =SIN(GAM  VHAD) 

CAvG=M/AR 

30  HEAD  {*>*201  ALPHAU*  jETA. FACTOR 
XFC  ALPHAU.GE.<»  ‘.  ) GO  TO  76 3 
IF  (FACTOR)  o ........  7 

FACTOR=0*4 
u 7 CONTINUE 
40  FORMAT  ' * 1 • ) 


50  FORMAT (/*14X* ’ALPHA  BETA  ’* 

+ * SWEEP’* 

1 7X  * ' AR  SPAN*) 

60  FORMAT (9X*6F10«2*F8,3*7X* II ) 

PI=3.14L592b5 

PGC0RR=1.0 

M=2*Nli-l 

NSSPAM-2 

BM=M 

MP2  = M+2 

MP3=M+3 

BMP1=BM+1 

BMP2=BM+2 

MUZER0=BMP2 

Ml -M+l 

NMMU=NiI-l 

CLAM-COS (SWEE P/RAD) 

SLAM=SIN(SWE'-  P/RAD) 
tlam=slam/clam 

63  DO  80  NU=1*M 

PH  I ( NU ) - { P I * NU ) / ( M+ 1 , ) 

F.TA(NU)=COS(PHI(NU)  > 

Y=ETA(NU) *8*0.5 
ABSY=ABS( Y) 

cau.  body ( alphas  * beta  » d i ab  » zw ing. bma jcr  * b - y , alfcou j 

C(NU)=GIRC(ABSY*YTAU,CTAU,NY» 1 ) 

I W-G IRC ( ABSY * YTAB * IwTAB*NY*l> 

ALPHA  (NU)=<  I W+ALPHAB+ ALFCOU) /RAO 
80  CONTI HUE 

ETA(MUZERO) =1,0 
ETA(BMP1 ) =-l .0 
PHI  ( MIJZEPO)  =U«; 

PHI (BMP1 ) =PI 
DO  41}  M=  1 * M 

84  DO  5u  NU=1 »M 

Bo  SUM=0,  I 

DO  20  MtJtO  P=1  »BMP2 

mu=mulo  p 

IF  (HU-1)  7 "5*8 ‘*105 

8 FBAR  ( N»  MlJ ) =0  , i 

90  DO  10  MUl = 1 . M 

FBAR  ( N,  MU ) =FBAR  ( H,  MtJ)  + ( MUl  *S  IN  ( MUl  ♦PH I (fj)  * C OS  ( MU  1 ♦ PH  I ( My ) 

1.  * ( 2 . / ( M+ 1 , ) ! 

10  CONTINUE 

105  IF  ( U—  1 ) 7/5*107,(3'' 

107  Tv.OAM=l  ,*<«*  PGC  OR  * T L AM* ( AUS ( E T A ( nu ) -ETA(MU)  /C(UU) 
T«0AP=1 . + { 8*  PGC OR-  '*TLAM*  ( ABS { ETA  ( HU ) • + E T A ( VU ) /C (NU! 

T rfOB-  PGC  OR  ‘ * i E T A { r iu ) -E  T A ( vy  j VC(UU) 

T«OC=B*POCOR  *(ETA(NU)+ETA(MU) ’ /C ( HU ) 

DEN0MP=1  . ♦ (U*  Ai»S(ETA(NU)  : ♦ETA(NU)  + PGC  OR  *TlAM)/C  ( HU) 
OEMOm'!=1  , * { ti*  ( AOS  ( ETA  ( NU)  -ETA (HU)  ' *PGCOW  *TLAM) /C ( NU ) 
F0UHA=: . ♦ (!j*A:sS(3TA(N>.i)  *PGCOP  *TLAM)  /C  ( HU) 

F OURH=  - * PGC  OR;j  * E T A ( f ,‘U ) / C ( HU ) 

FlVEP=(SQRTf  'Ta'OAP*  2)*(TaOU*«Z)  /UENOMPI-I.u 
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FIVEM=  ( SQRT  ( ( TWOAM*  • 2 ) + ( TWOB*  t-2  ) /DENOMM ) -.1 . U 
SIX=2.0*TLAM+SQRT ( ( FOURA*  2 ) + { FOURB* *2 ) > / < DENOM^*OENOMP) 

IF  { ETA  { MU)  1 0,12'), UO 
110  IF(MU-NU)  1 5*1  8 . 1 8 

115  THRE  =1 ./ ( B* ( ETA ( NU)-ETA ( MU)  ) -PGCORi  VC ( NU)  : 

HSTAR(NU»MU)=THRE  *fivep+six 
GO  TO  130 

118  HSTAK(NU.NU):-TLAM+SIX 
GO  TO  130 

120  IF(MU-N'J)  125  * 128  * ' 25 

125  THRE::  = 1./(B*(ETA(NU5-ETA<MU}  j*PGCORR/C(NU) 
HSTAR(NU»MU)=THREF*FIVEM 
GO  TO  150 

128  HSTAR(NU#MU)=TLAM 

ETA  BAR:  ETA  .GE.  ZERO 

130  CONTINUE 

IF (MU-BM) 180*10 ) *20 

180  SUM: { -1 , / ( 2 . * { BM+ 1 ) ) * { FBAR ( N. MU 5 *HSTAR ( NU* MU)  ) +SUM 
20:  CONTINUE 

GBAR(NU»N):(-1./(2.*(BM+1) ) )*( ( FUAR ( N* MUZERO) *HST  AR ( MU * MUZERO ) 
+ +FBAR(N*UMP1)*HSTAR(NU*0MP1)  /2.01+SUM 
30  CONTINUE 
40  CONTINUE 

U0  60  NU:l*M 
00  405  N: 1 * M 
LB ( NU*  N ) :0 . I 
DO  402  MUl - 1 »M 

LB ( NU f !) :LB(MU# N) +MU1 *S IN ( MUl * PH I ( NU) ) *5IN( MUl *PHI ( N) ) / < M+l . ) 
1 *5IN(PMI(NU)  * (2 ♦ ) 

4 D2  CONTINUE 
4)5  CONTINUE 

DO  50  N : 1 • m 

4 30  A(NU>N) :LB(NU*M)+GBAR(NU»N) *U * P&CQR'  VC ( NU ) 

50  ■ CONTINUE 

A ( I |U*  N"1  :A  ( IJU*  MU) -H/  ( C ( NU)  *C'_Am*  PI ) 
nO  CONTINUE 

IF  (PRINT)  7P  * ' 'HO*  ' >7  . 

«-,RITF(6»  ( NU*  N.  A ( HU* ! I)  > U=1  * m)  . NU=1 . M) 

oO  FORMAT  ( 4 • . II  . » ) :' .FI 'Jv  5) 

W R I T E ( 6 * 973)  ( N»MU,FBAR(M.MUi  . N:\.M)  .mij=1,mP2) 

HU  I T F ( o » • A)74  J ( NU*  M! ! . HSTAP  ! f JU . V !)  . NU:  l ♦ M)  . MU=  1 * MP2  ) 
n'?73  F0«mAT(4(*  FBAR(  • » II  . * . • » II  » » ):•  * F1C.5)  ) 
t>'»74  FORMAT  (4  ( • HSTAR(  * * U ■ ' > ' » II  * ' ) = * * F>  0V5  ) 

WHlTE(t»*  -970)  ( NU*ii.  LbdlU.t,')  *lt')=l  *M,  .N- l*M) 
n«70  FORMAT  ( * LlH  * » 1 1 * • ' * 1 1 * * 1 : ' . FM.5 ) 

WRITE  (6*6972  ) ( ( NU*  M.  GUAR  ( NU.  N)  ,nu:i*M'  »*t=l«v) 

- >*97  1 FORMAT  (4  ;•  u -i  * . II  * • . • . It  * * ):» .Fn.5’ • 

6072  FORMAT (4 (*  GUAP(»»I1**  » » 1 1 * * ) :*  *F8.5 ) 

7**«0  CONTINUE 

eOK=AHS(EOKPCT/lO  ,#ALPHAH) 

U C £0^  . IT  . .0  5)  f.OK=.0‘»5 
LINE-I 


WRITE (6*7982 ) ALPHAB 
7982  FORMAT! *1  ALPHA: FIO ,2 ) 

CAL  GUES;(N  :» FACTOR*  LINE) 

CALI  SOLVE <NU»EOK»  IT  1 ER»ESTAR» FACTOR) 

WRITE (6 1 910)  IT TER » : OK»  EST AR 

910  F0RMAT(1X# * IT  ER=* » 13# * EOK= * , FI 0 . 5 » • ESTAR:*  »FlU«b) 

DO  020  1:1  * M 
DALFE:ALPHAE( I)*RAD 
CL1=2.*BxG( I)/C( I) 

DALFA=ALPHA  •( I ) *RAD 

920  IF ( PR  I NT • GT • 0 . 0 ) WRITE(b»92b)  ETA { I > , DAlFE * CL1 * UALFA 
925  FORMAT(lX»4F10.b) 

WRITE (b  »4U ) 

WRTTE(6*1  Is)  (TITLE{K  K)#K  :=1»  8) 

11.1  FORMAT  (/tlX*  8A10  ) 

11  2 FORMAT (8A10) 

WRITE (6.50 ) 

WRlTE(b.oO)  ALPHABf  >ETA»SWELP» AR#B 
10b)  CONTINUE 

SUM48T=0.  i 
DO  120')  NU=1#M 
SUM48(NU)=0.  1 
DO  1180  N=1*M 

SUM48 ( MU) =SUM48 ( NU ' + LB  ( NU » N ) * G ( N ) /2  • 0 
1180  CONTINUE 

SUM48T=SUM48T+G(NU) *SUM48( NU) *SlN(PHI (NU) 

120  CONTINUE 

CDI921  = (PI*L!/CAVG)/(M+1  )*SUM48T 
1250  DO  U'J  1 = 1. M 
MP2MI  = MPi-I 
I PI  — I 4 1 

PHI ! ( ! PI ) =PHI ( I) 

PHI  12  (I) = PH I ( 1 ) 

ETAI!  { IPl)=-ETA( I) 

Y = £TA(  I)*B‘.b 
AlJSY  = AfiS(  Y) 

CMU=GIWC(ABSY»YTAB,CMUTAB.NY.l) 

AE= AE/CLAm 
CAMAOVt I P 1 > = G ( l:*U 
r 0o=o . 

ALPHA |11l =5UM4B ( I ) 

CL~(2.!)*B*G(  M /cm 

EPS:? IRC ( Y* YTA8  *EPSTAB.NY« 1 ) 

HLAM:v.0S(  VwE  P-LAMHLi  /RAO)  /COStLAVHL/HAD) 

0ELFS=G1!,C  { Y»  YTAB.DElF  TO.  NY * I 1 
0ELF=ATAM(  'AritUELFS/RAOJ  •RL,AV)  *PAD 
OELJS=GIHL{  t » fTAR.uM-LJTb.NT.  1 ) 

OF l.  J: A T A r j ( f AN;  ■ 0ElFS*0EL  JS) /WAD?  *HLAM)  .RAU-OtLF 

CUPOLN=l  *0-«  SIM  ‘ EPS-LAMML)/RAU)+C05(EPVWAD)  IN( LAMHL/pAD) 

1 * ( l ,0-C0S(  ;OLLF*'-'ELJ)/“ALM  *1 

C MUP= C U P OC  U ♦ C i 

CLi  0 AUPl)  = (CL*CNajP#SI»K  :ka'AkE)*(  ALPHA  f I ) -ALPM ( 1 ) #C  ( I ) 
1 / CAVG 

CKOL  ( I PI  1 = (-CL*C  ! D/CAVG)  *V/B 

U9 
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CDICOC  IPl)=CLCOCA{ IP1 ) *SIN ( ALPHA A < I ) -ALPHAE ( I ) 

CM-O. ) 

ALFAE1 ( I)=ALPHAE( I) 

ALFA  I ( I)=ALFAI1  { I) 

130  CONTIMJE 

CLCO'J!  (MP2)=0.0 
CROLL 1 MP2 ) -0  « 0 
CDiCO'  (MP2)=0. 

GAMAO /( MP2)=0.0 
F * II 1 ( ^tP2)=PI 
EVA I (MP2)=1.0 
CL Cr  ,A{1)—  0.0 
CKLL.(1)=0*-J 
CL  ICOC(1)=O.U 
GAMAOV(l)=u 
ETA  I [ ( 1 ) --  .1  • 0 
PHI  I i 1 )=0. 

WRITE(b#7/Y  ')  <Kk  »PHII  IK'  ) .K  K=1*BMP2) 

WRITE  (o*7  8 ) ( KKk  * ETA  I I ( K>'  ■< ) * KF  rl,BMP2) 

WRI TE(6*7 7 ’9)  (K--*CLCO  A'k  ■ )*K  • =1»BMP2) 

77  F ORMAT ( 4 ( * PHI i ( * » 12* • )=’ *FlU.b)  ) 

7/78  F0RMAT(4(*  ETAI  I ( * * 12* ♦ ) = * *F10,b) 

7/  9 FORMAT ( 4 ( * CLCOCA (♦* 12 *•)-•* F1U , b ) 

CMUT=0.  l 

DO  130b  KU=2»NY 

KUMl=KU-i 

CMUT=CMUT+ (CMUTAUI KUM1 ) *CTAB(KUMl ) +CMUTAU ( KU ) *C TAB ( KU > ♦ 
1 (YTAB(KU)-YTAB(KUMl ) )/ (B*CAVG) 

1305  CONTINUE 

CLT=QTAB(CTAI I *CL  0 A*1 *MP2 )/ '.U 
CROL  T=OTAB(ETAI : *CROL'  *1 ■ MP2 ) / ’.0 
CD  IT  = QTAii  ( ET  A I i , CD  I - OC » 1 * MP2 ) / .0 
SUM4o=U.U 
DO  1310  N- 1 . M 

SUM4«j  = G(M)  ♦SlfKPHI  CN)  +SUM4« 

13 '0  CONTINUE 

CL921=(  !PI*B)/(CAVG*  (M+l  ) *{SUM4e>) 

C 

C CX  ALClILATION  RER  L.  BARNLTt 
LGAMA=U.921  **2/  <PI«ah*CI)I»21  ) 

AU4.*CL921/(PJ*CLC0  A(f|.) 

XKHAKT4  « / ( P I * 2) 
h 1=  ,U-(EGAMA*XKUAK/ (u , *A1 ♦ • 2 ) 

£i  2 = 1 .0<M2.*ETACMU*CMIJT*XK»AP*G  1/(PI*AP*A1*  /)> 

O 3-CL T * 2/(PI*AR*EGAM  *0  2) 

CXT  = ET  At.MU*C?*«lT-CD0T-G-3 
C 

WRITE (b*  710  ) 

710  FORMAT  (/*1X.*  CL  CL'CIWC]  CDP** 

1 • CX  ETA  CMU  CKN<OL«_)*’ 

720  WRITE  U>*73**»  CLT*  L92  l*CUQT*CXT*ET  ACMU*  C’>AJT  • CROl  T 
730  F0RMATUX*4Fl0.5*Fft.3*2X*.:F]0.b) 

740  t.N ! TE  (o * 7b0  ) 
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750  FORMAT  {/»bX**  ETA  GAMA/V  CL+C/CAVG'* 

2 * ALPHA  E • ) 

OUT=0.  ) 

ABSiiET-AiiS  ( BETA ) 

DO  7 2 1=1*20  i 

YW=-1.0+DELY*< 1-1 5 

IF{  ADS  '.ET.LE.  Oil  ) YW= J . O+DELY*  ( 1-1 ) 

CHECK=1.-YW 

IF  (CHECK)  7 75  * 756  * 7 j 
756  0UT=1 • 0 
YW=1 .0 
YViPHI  = 0 • ) 

GO  TO  1759 

757  IF  (YW)  759*768.759 

758  YvVPHI  = P 1/2  • 0 
GO  TO  1759 

759  IF ( YvV+1 ) 775*1756*1768 
1756  YrtPHlrPI 

GO  TO  1759 

1758  YvVPHI=ACOS(Yw) 

1759  GAMAW  = G IRC ( YwPHI * PH  I i * GAMAOV  * MP2 . ) 

C ■•,‘RITE  (6*925)  * Y,V*  Y.nPHI 

CLW=GIRC  ( Yv/PHI  • PHI  I * CLCOCA » MP2  * 2 ) 

CDI  //=GIRC  ( Y »VPHI  * PHI  I . CD  ICO  - * MP2  *2  ) 

ALPW  = GIHC ( YWPHI * PHI !2 . ALPHA  I , M/2 ) *RAD 
CM//=G IRC  ( Y //PHI » PHI  I * CMC 0 * MP2  * 2 ) 

ALPE=6IRC(Y.vPHI/PHI  ! 2 * AlFAEI  . m»  2 ) *RAD 
TEMP5  = GIRC(Y.VPh;.PHT  12*  ALFAI,M  2)*RAD 
7o0  v '/RITE  (0.7  0)  YvV.GAM-w/.CL.v  * ALPE 

770  FORMAT  ( 5X » 1 F 1 0 . 5 » iX * IF  0 .5 . 5X  < L f 1 0 . 5 . 5X  * IF  • u . 5 . . 

1 FLU. 5) 

IF  (OUT)  7 '2#  '/2.  r 5 
772  CONTINUE 

7 8 IF(flPTS)  782*781.781 

781  CAL  FIELDfNN.riPTS.B.ALPHAH.  (ETA.SrtE  P*  GAM  > 

1 .PHI  I*  GAMAQV  *XP‘J*  YPD*2PU»MP2) 

7 ft  2 I F ( N r-  i 6 ) 5 li » i > * 7 6 } 

78J  STOP 
END 
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SUBROUTINE  GUESS  ( NS' >*  FACTOR*  LINE) 

COMMON  ;)SMPAN*ti*  ALPHAA<5  )*A(id'>  )*C(5  )* 

SKIP  (5  ) *X  (5b  *5  ■)  * ERROR (35*  .j)*E(3  )*XSET(5  ) 

* SK IP2 (35) * YTAH (2U ) *ALPZL(20) 

* NY* CLAM 

DIMENSION  TEMP(5  * i ) * G ( 5 )*ER'SET{5  ) * TEMALF  ( 5 i ) 

RAD=b7.296 

PI=3.14!b926b 

NCP=NS 

M-2+NS  -1 

IF (NS  PAN. EG. 2)  NCP=M 
NCP  1-MCP+l 
IF  (LINE-2)  40  * 2 -b  * U 

TN IS  PORTION  COMPUTES  A LINEAR  SOLUTION  FOR  A FIRST  CUES  .. 

00  lb  Cl  1=1*  NCP 
Y = AJS  ( (*.b*COS<PI*I/(M+l) 
alflu=g  I PC  ( Y * Y TAB  * ALPZL*  MY  * 1 ) 

TEMALF  ( T)  = ALPHA  '•  ( D-ALFl.O/RAD 

do  iom  j=i *ncp 

IF(I-J)  bU  *b  1 * bb 
TEMP ( I, J) = A{ I.  J) 

GO  TO  10  i 

temp ( i*  i } =a ( I*  n+d/(C:  I)*pi*clam) 

CONTINUE 

CONTINUE 

CAL'  GJRV  ( TEMP*  NCR*  5 .0.')  1*1  EW  ) 

IF  ( IIP  ) lbU  • 17  1*1/  • 

aR  I TE  ( o * 1 ’)b  ) IER 

FORMAT ( IX* • IER  = * * 12  * • IN  INVERSION  FOR  LINEAR  SOLUTION*) 
CAL  YATvpy  ( TEMP,  TEMALF  *G.  NCR*  ICR»1*5  * » > ' ) 

00  20  1=1. NCP 

y C I * 1 )=rf»G{  I)/(C(  I)  *PI«CLAM)  t-ALF!  U/RAU 
XSETt  I)=X(  1*1 ] 

■vR  I TF  t o * l MO  ) I.X(I»1) 

FORMAT ( i X . * 1=  * . 12 . • X ( 1*1 ) = * * F l U . b ) 

EP  0 M 1 * 11  = 1 .0 
xnCp=n;  •’ 

kp  o -d.NCP  n=i  .u 

CAL'  CALC (NS  *ER  SET*  (1) : 

xsuM=o.t-‘ 

DO  2 u 1=2*  NCR  1 
1 vl =1-1 

EP-Q  - ! 1*1!  = t>  - SET  5 I d ) 

continue 


THJS  PORTION  MAKES  NCP  MORE  GULS‘  L B&SLO  ON  FIRST  GOES  * 

00  50  0=2 » NCP  1 

00  2Sb  1=1 ♦ NCP 

IF  ( J-NCP  l I 244*25  s»25'> 

I F ( 1 — 11  25b  * .) 5b.  74  0 
X( I.u)=FACTOP*X(  1*11 
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GO  TO  248 

240  IF(J-I)  245  > 235 : 245 
245  X(I» J)=X(I»1) 

248  XSET ( I ) -X ( I » J) 


250  CONTINUE 

CALL  CALC  ( NS',»  ERRSET  » E(  J)  ) 
DO  260  1=2 1 NCPP1 
I Ml- 1-1 

ERROR  ( I » J)  =ER  'SET  ( IMl ) 

260  CONTINUE 
300  CONTINUE 
50  *.  RETURN 
END 
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SUBROUTINE  BODY ( ALPHAB# BETA , DIAB, ZW ING, BMAJOR* B , Y - ALFCOB) 
IF(DIAB)  50*5  ' >25 

25  IF(BMAJOR)  30,50 >26 

26  IF ( BMAJOR-DIAB)  60,30,60 

CIRCULAR  BODY  CALCULATIONS 

30  TEST=Y*fc2+ZWlNG**2-(DIAB/2.0)*  2 
IF  ( TEST  ) 50,3  >,55 
35  ETABOD=2.0*Y/B 

ANUM=(ETA30D*B/DIAB)*  2~(ZWING*2./DIAB)*‘2 
ADEN=(  ETAB0D*B/DIAB)*+2-(ZWING*2./0IAB)*  <2)  i=  i 
+ +4 » 0* ( ETABOD* ( B/DIA J) * ( ZW ING*2  « 0/DIAB ) ) * -2 

BNUM=8.0*ZWING*Y/(DIAB*  2) 

BDEN=(4.0*( (ZWING/DIAB)**2+(Y/DIAB) *f2)  ' 2 
ALFCOBr ( ANUM/ ADEN) *ALPHAB+ ( BNUM/BDEN) *BETA 
GO  TO  10 

NO  BODY,  RETURN  WITH  ZERO  ALPHA  EF  ECT. 

50  ALFCOB=0.0 
GO  TO  10 

el' iptic  body  calculations 

BMAJ0R  (VERTICAL)*  DIAB  (HORIZONTAL) 

60  Al-DIAB/2.0 
B1 -BMA JOR/2 • 0 
A2=(  ( Bl/A.L  ) * *.2-1 . 0 ) /4 . 0 
R2-  ( Bl/Al)  + l«(J)/2»U)*t2 
TETA~ZW  ING/A1 
TPSI=Y/A1 

TEST-  ( Y/ A1  )*  i-2+  ( ZW  ING/B1  )*  *-2-1.0 
IF (TEST)  o2,o3*o3 
o2  ALFCOB=0.  > 

GO  TO  10 
63  I T -0 
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648  DO  75  1-1 » 80 
IT=IT+l 
ET  A-TET  A 
PSI=TPSI 

TPSI-Y/A1+ ( (A2*PSI)/(  PS  I ) * *2+ ( ETA ) + » 2 ) i 
TETA=2WIM6/Al-( ( A2*ETA)/ (PSI**2+ETA*  2 > 

Ic(  ABS(ETA-TETA)-.0  "'1)  6b»  » )»7-i 
65  IF  < ABS(PSI-TPSI)-.00  1!  8)*80>75 
75  CONTINUE 
80  CONTINUE 

C WRITE (6 » >9°)  IT»  Y > ETA » TE  FA>  PSI > TPSI 
ETAzTETA 
PSirTPSI 

C69  1 FORMAT (1X» I4.6F10.5) 

BNUMz  { R2  *■  A2  ) * J « *PSI*ETA 

80EN-  ( PSI*  2-(ETA*  2M-A.')*  .’+ ( 4 . * { PSI*  * 2 ) M E i‘A*  ♦ 2 )) 
ANl)Mlz(PSI*ti>+ETA*!2)*!2+(PVI*t  2-ET  At  : ) ♦ ( R2  + A2 )*R2* A 
ADEU1=(PSI*  2-ETA*  r2  + A )*  2+4.MPSI*  2) ' (ETA*  2) 
ALFCO«=  ( ! ANUM1/A0EN1  )-i  ,U!  *ALPHnii+  ( UNUM/LJQEN)  +UETA 
10  ! RETURN 
EMU 
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SUBROUTINE  GJ^Vl  AfM  . NL#  EPS  EL  » ItRR ) 

A IE  INPUT  AR*AY  WHICH  WILL  CP  DESTROYED 

n is  dimension  cr  a 

CPC  I L 1C  ICC  T VALO'r  FCP  eiVfT 

IERR  IS  VALUP  RETURNED  IF  PIVOT  IE  TCC  EMAIL 

DIMENSION  A (1) #3 (125) »C (ID 6) tlFCiDS) »IQ(123) 

ICRR-0 

DC  140  K z 1 » N 
«IVCT=0.0 


120  T=X*N 
DC  2 JzK.N 
IDFXz  ( J-l ) *NL  *1 

IFC A£S(A(ICEX})-A:D(FIVC’ ) ) DtZf : 

C^NTINU" 

pIVCT  = A (HEX) 

Ip ( K ) zl 
IQCIOrJ 
CONTINUE 
CONTINUE 

TF{  AOSr’IVuD-EPClDnCtlUC'.: 
CONTINUE 

Ic(  IFUO-K  ) 4 1 C » 4 

CONTINUE 

SWAP  ROWE 

CO  i J-  1 • N 

InX ZIP (K  ) 

IDFX-  (J-1)*NL-Hr  X 
KCFXZ ( J-l) • NL ♦ K 
Iza(IDEX) 

A ( IDEX ) z AIKCCX  ) 

A( KDEX ) zZ 

CONTINUE 

CONTINUE 

:<-(  IC(KI-M  )?tTf7 

CONTINUE 

swap  column: 

CO  9 I z 1 * N 
i^yziDr/ ) 
icrxzd^x-i 
KCCXz  (K-l ) «nl  * i 
EzAtICTX) 

A( TDCX  JzAtKCEX) 

AIKDEX) zZ 
CONTINUE 
CONTINUE 
OC  12  JZl.U 
KDrXz|J-U«M*K 
JDEXzIK-; *«NL*U 
If  t J-K ) lit  10  1 11 
CONTI NU“ 

01 Ulzl. /PIVOT 
0 ( J ) z 1 . 

OC  TC  1- 


11  CONTINUE 
S(Jl=-AtKOEXJ/riVOT 

C( J)=A( JDCX) 

12  :CNTlNUr 
AtKDEXI  -O.C: 

At JDEX J^O.C 

i ? CONTINUE 

r*  14  1-ltN 
DC  14  J - 1 1 N 

idex-u-i ) 

A f TCEX  J - A ( IC  EX) ♦C(I)*2  ( J) 

14  continue 

14 l CONTINUE 

DC  DU  KP=1»N 
KrN* 1-Kr 

Tc( IP(K)-K)15»17#15 

15  CONTINUE 

DC  1?  I-liN 
▼ p y r I T { X } 
TCEX-(InX-l)*NL*I 
Kcrxr  IK-1 
Z rA  t ID  EX  J 
AtlCEX 1 -A  t K C r v > 

AtKDEXI-I 

ic  continue 

1?  CONTINU" 

:rcccx)-KJi?»ZL»it 
19  DDNTintr 

DC  IS  J=1«N 

:px=:ot*  ) 

:;ry;(J-D**-L  + lf'X 
NO" X - l J-  1 J • ?»’_ ♦ K 
::Micrx) 
a?td:x)-aikdex> 

AtKCfXl=: 

13  CONTI  Nt!" 

:c  cc»:t:nu': 

cc  t:  :i 

l(.t  CONTINUE 

:-0tDr-l 

;i  CONTINUE 

p-TTU^N 

tv: 
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FUNCTION  CIBCURC.X  .Y.NX.IC) 

AITKIN  INTERPOLATION 
AP  C = IN?rD£N PENT  A°f UME NT 
X ^INDEPENDENT  TAEL! 

Y=CErEN3ENT  TAELE 

NX "NUMBER  VALVE3  X TABLE 

IC  = 1 FIBCT  CPDER  IN TCPF CLAT TCN 

10=2  2E00NC  o9zzr  interpolation 

2I*/  ENEICN  X ( 1 5 » Y ( 1 J *XX (4) * Y Y (4  ) * Z E ( 2 J*FFC2  J 

TEN0=NX-IC-H 

I F ( X C 1 J - X ( 2 ) Jir.2r.2L 

A3CEN0IN3  CREEP. 


ir 

2C  15  1 = 1.11 
IF  ( X t IJ-A.BG 

15 

CONTINUE 

IE 

:_,=nx-ic 

IFXrO 
30  TC  4B 

2C 

DC  25  1 = 1.1! 
IF(  XtD-ARC 

2 r 

CONTINUE 
00  TO  lc 

:r<:c-i)25» 

7 - 

I'C 1-1)23,3 

3E 

!■»  = ! 

00  TO  43 

27 

T — T _ *-  £ 

TO  43 

4C 

:ri:-2)4m 

41 

I"X=0 
GO  TC  3r 

42 

:*"y  =i 

Cr  T r T ■* 

45 

^0  50  1=1.4 

xx  tr  >=x  c:  i 
YY(I);YC3> 

5C 

:r=:r*i 

30  GO  1=1.3 

cc 

EDI)  =XX  (:■» 

*1  « Tf> 

7r 

rr  t : ) =r r t ; ) 
03  30  :=i*2 

:r « i » = < yy  1 2 
:r(  :c-mou 

Ef 

CONTINUE 
3*i  : = *.: 

?t 

:n:i  = mii 
:ri  :ex)ioo» 

or 

* v 

Ercnnrm 

IOC 

oipo=:ri d 

•JETUPN 

~ \ ** 

. % . 

c 


SU3RCUTINZ  KA’KPYtPatfPCSTfPRCS*  T,U»K»  M?R»MPC»MPO) 
DIMENSION  PF.E  « 1 5 * PCSMD.  PRCO(l) 

3CU3LE  "RECISEON  3UMMY 

CRT "RE-MUt.TlFLlEF 

'’C3T--  PCBT-M'JLTIPLIER 
P°CD—  PRCLUCT  MATRIX 

I --  NUMBER  oF  ROWS  IN  P">.r  t NUMBER  OR  ROWS  IN  PROS) 

J — NUMBER  CF  CCLS  IN  Pr  t ( NU  M3  E R CF  ROWS  IN  PCST) 

K --  NUMBER  CF  EC  LS  IN  PCBT  (NUMBER  CF  CCLS  IN  pOST) 
MrR  --  ROW  CIVENEICN  Cr  &•:£  (FIRE"'  CUBE C RT pT  ) 

MPC  --  ROW  DIM  EN  3 ION  CF  «’3ET  (FIRE!  -»U°S:RIPT) 

MFC  --  PCW  DIMENSION  OF  PRCE  (FIRST  SUBSCRIPT) 

30  10  L - 1 # K 
UP  AP  TrM-3  *L-M"S 
LPCRT ~y.nZ*  (L-I  1 


COMPUTE  PRODUCT  EY  COLUMNS  rp CM  SUMMATION  CF  PRODUCTS 
C C F RCW-"RE  .'IN"  OCLLVN-  p'>"3  ELEMENTS. 

30  m M=1»I 
U°ART-L'>C.RT 


If  AR  T z* 

LPART -Ln ART ♦ ! 
"UMMY-C .030 
3C  3 N- 1 # U 


JrART-UPART-»  1 

DUMMY  rP^E  ( IP  ART  )*PCST(  UP  ART!  ♦DUMMY 

IPART  = I,iART*MpR 
PRC  D (L°  ART ) -DUMMY 
RETURN 


R 

1C 


C CT  A 2 IN -2  CP  AT  ’CM  PCUTIHE. 

C CALL  INC  StQUENCr 
C ANC-CTAC(X#rX»I*j) 

C X i:  A TA2L-  ZF  V VJES  OF  T HE  A^CjMENTt  FX  13  THE  T A2LE  OF  VALUES 
C CF  THC  FUNCTION,  rx  III  - r(X(in.  THE  IMLC^ATICN  IS  F ROM  XII) 

C TO  X( Z) . 

C A F AH  ABC  LA  IS  FIT’ES  THE  POINTS  (EY  THREES  ).  ^HZ  AREAS  RESULTlNS 
C ARE  SUMMED • 

r 

rUMCTIC'«3T  A?(X»FX*I»J) 

DIME N SIC NX(1)  try  (1 J 

:um=o. 

M- J-I-l 
IF { M J 5 1 4 # 2 

n i — |—  p 

DCTKrltLtE 
TI:FX  (H-tl  J-FX(K) 

T2 -X («♦! )- X ( K ) 

T3  = X CK > -X  C K J 
T4=XtK*?)-XU*l) 

AI=T1/T? 

A2  = (FX(K»2)-rv(K)-AI*T3)/rp3*T4) 

C2  = « Al-  A2*  (X  ('■')  *X  (Xn))  1/2  . 

C?rFXlK)-Al*‘X(K)*A2*X(K)*X(X«-l) 

XT2  = X tK  )*«2 
X22 -X  C K>2  ) **2 

xpi  = xo:*x (k  i 
X23-X22*X(K*21 

: :jv=:uy*{AL*  (y:3-yr2i/7.*c:*cxcr-xt2i*c?«t»(K*2  j-xiKDi 

T F t MCS (M«  2 ) ) I#4»E 

4 SL’M=:UM*t  <X{  JX-XC  J-ll  J*  «rX  < J)  ♦FXlv-1  J )/2.  ) 

5 T’A2=SUM 
TETUFt; 


SUBROUTINE  SOLVE ( NS' .» EOK» IT i ER»  ESTAR*  FACTOR) 

COMMON  NS- -.PAN*  d »SKIP(1  8U»X(3  > 33 ) t ER  'OP  (3  >3  ) » E ( > i ) * XSTAR  ( 3 ) 

DIMENSION  ERi'.STR(3  >)  »B<3  )»  i )»  UNITY  ( 3 .*  ' } » DXSTAR { 3 ) » DER-GT  ( 3 >) 

ITTER=1 

NTRY=1 

NCP-NS  : 

M=2*NS'.-1 

IF  (NS  .PAn.EQ.2)  NCP-M 
NCPIJ1  -NS  P+ 1 
RAD=b7.296 
b . CONTINUE 

DO  lbU  J- 1 » NCPr’l 
DO  10  i K=1*NCP->1 
B( J#K)=ER  0 ( J#K) 

1)  ! CONTINUE 
150  CONTINUE 

CALL  6JRV ( B * NCR  1 ♦ 3 -.0  1 1 IER  ) 

IF  ( IER  ) 160  » 170  > 16') 

160  WKITE(b»lub)  IER' 

16b  FORMATdXf  ♦ IER  I4» ’FROM  GJRV  MATRIX  INVERSION*) 

IF (IT  ER-10  ) 9bDrl64»lo4 
164  DO  nO  1=1# NCP 

qO  WRITE  (6#  1 ) (X(  I#  J)  # J=l»NCPuU 
bl  FORMAT (IX*  *X( I * J) * » 1UE1  .4) 

•<;RITE(n*  2)  (E!  I)  * I=1»NCP  1 ) 

62  FORMATdX*  ’E(  I)  • r 10E  .4) 

DO  70  I=1»HCPN1 

7«)  «RITE(o*71)  f ER  O'd  I*  J)  #J=1  *NCPJ1  ) 

71  FORMAT ( 1 X # ’ ER'OR ( I » J)  * » 1 OE  ;.4) 

IAL  MATMPY(ER  OR»  B » UNIT  Y » NCR  '1  * NCPm  l # NCP-  l #3  » '•  #3  ) 

00  In  1 = 1 » NCP- ‘l 

16  '/R  I TE  (o  » 1 .7  ) ( UN  I Tv  ( I * J)  * J=  1 * NCP  1 ) 

1«7  FORMAT  (lx*  • Nfj* TY  ( I • J)  * » IUE  .4) 

DO  loH  1 = 1 »NCP  '1 

InH  ,\R  IT  E to  * l .9 ) (U  I#  J)  # J=1  *flCP  1 ) 
led  FORMAT ( IX* »BINV( I* J) * » lUEi  .4) 

GO  TO  950 

l?u  DO  y<'-  I * = 1 * NCR 

SUM= 

('0  20  U.dfHCR'U 
SUMrSUM+xt I ! » J..U  *B( j » 1 5 
20  CONTINUE 

XSTAR  (I  • ) rviM 
0XSTAr(  I I )= XSTAR ( I ! ) *RAQ 
ji  CONTINUE 

32U  FORMAT? 1 x# »• STAR* » 10E ; .41 

CAL  CAL,  (NS  *VR -STR.ESTAR) 

DO  325  1 = 1 * NCR 
OtR  ST(I)=ER  STR(  n*RAU 
325  CONTINUE 

3 4 FQPVAK1  <»  ’ER  ST’  ♦ lOEl  .4) 
rtR  I TE  ( <> » 30 1 ) 65  TAR 
39n  FORMAT  I IX#  »E STARS’  #£14.7./  ) 
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591  FORMA  (1X»*ESTAR=»,E14.7) 

IF(ESTAR-EOK)  400#4Jj»bi.) 

40o  GO  TO  10  » i 
5JH  EMAX-E  ( 1 ) 

NEXCH=1 

DO  70  K=2#NCPU 
IF  ( EMAX-E ( K)  ) 6Uo » 7 ,»7-Jl 
60  NEXCH=K 

EMAX-E ( K ) 

70  * i CONTINUE 
EMIN=E(1 ) 
nmin=i 

DO  726  K ~2  t NCPi’l 
IF  (EMIN-E(K)  72b»72b*71> 

716  EMIN=F(K) 

MMIf.JzK 

72b  CONTINUE 

IF (EMAX-ESTAR)  7b0»760,80 
7b0  wRITE(b»7b j)  EST  AP , EMAx 

7b  FORMAT ( 1 X » * ESTAK- ’ » F8 • 4 # * IS  GREATER  THAfl  EMAX= * > Fti . 4 ) 
GO  TO  960 

760  WRITE (6 » 76b)  EMAX, 

76b  FORMAT ( 1 X » * ESTAR  IS  EQUAL  TO  EMAX= ♦ » Ftt.4 ) 

GO  TO  9bU 

HO  DO  HbU  I 1-1 » MCP 

X( I ! *NEXCH)=XSTAR( I ) 

HbO  CONTINUE 

DO  60  JJ=2fNCP.'l 
J .'Ml- J - l 

ER  OR( JJ»NEXCH)=ERfSTR{ JJM1  ) 

90  CONTINUE 

E ( NEXCH) :ESTAR 
IT  L 6: I T ! EP+ 1 
GO  TO  bO 

9b  0 WRITE  (»»9..0  ) E(NMIM)#E0K 

9o0  FORMAT  ( - v I N RMS  ER  0 =*  *F'lU.b»  • LOK=*»FlU.b) 

9o2  Kr  =2 

IF  (fJTRY.EQ.2)  K*  =1 
IF  ( K ) 962 » '*<jJ  » 6t>H 
‘M7  FORMAT  III!) 

96  H GO  TO  ( 9HU»9o>*  U 1 1 ) #962*  ’>2#  * »s  ) » K>  -- 

C 

C K<  ■ =2  at  .IV  CAL  CUES.  A&AJN  USING  THE  . M J n CASE 

C 

9u  DO  » 7 1 I ; -1 , MCP 
I PUI  *1 

xU  ♦ i ) -x  ( I ; * mm  is) 

ER  O'  { Is  PI  » > -ER  O-OI  IPltMMls) 

971  CONTI HUE 

E ( 1 j-Einvjm) 

FAC  T OR-O .ft 
MTRY-2 

CAL.  GuES  IMS  ..FALTORU) 

GO  TO  SO 
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C 

C KKK=1  w'E  want  TO  PRINT  OUT  MIN  ER:'OR  solution. 

c 

980  DO  982  1 1=1* NCP 

XSTAR ( I I )=X( II»NMINJ 
982  CONTINUE 

CAL!  CALC  ( NSG* ERRSTR * ESTAR) 

GO  TO  10 
C 

C KsK=i  PRINTOUT  MIN  X AND  ER^OK  SOLUTION. 

C 

9 >0  DO  9 4 I i =1  * NCP 
HPi=ir+i 

DXSTAKI I i )=X( 1 1 » NM IN) *RAD 
DERUST (II) =ER  - yn  I j PI  * NMIM) *RAD 
994  CONTINUE 

WRITE (6* 320 ) ( DXSTAR ( I) . Irl *NCP) 

wR I TE (b  * A 0)  (DER  ST(  I ) # 1 = 1 ,NCP) 

WRITE ( 6*  390 ) E ( NMIM) 

GO  TO  9o2 

K < =4  CHANGE  ONE  OP  THE  MIN  X VALUES. 

I O'  WRITE  ( b ► 1 ! ; 0 ) 

1 U FORMAT!  « TYPE  <I2*F10.  ) MO.  AND  ME/,  VALUE  FOR  X.») 
REA0(b#l  b)  lN»XNEw 
IF  (IN-NCP)  1 7*1  7*:  0 

1 5 FORMAT ( 12* F10.  ) 

1 7 X(  IN*UM1N)=XNEW/RAL> 

DO  1 20  I . =1 *NCP 
XS  I AR  ( j J ) =x  ( I * r im IN) 

l 2u  (.out  i hue 

CAL  CAL  ( NS‘>»  ER  STR,  ESTAR) 

E ( NMIN) =ESTAm 
DO  I 2b  I i = 1 » NCP 
I P1=I •+! 

ER  OR  ( I ! PI  * NM  I N)  = t R ST  -’  ( I ! ) 

OXSTAR  ( I J ) = XSTA*'  ( I )*RA(J 
UER'  STf  I I )=FR-  ST  '(  I )*RAQ 
I 2b  CONTINUE 

WRITE  (t,.02u)  ( OXSTAR  t I)  » 1 = 1 . NLP) 
wR  I T (n  .*  A si)  (OPR  ST  ( I ) . 1 = 1 . NCR? 
wRITE(b*A9iJ)  E(MMI»l) 

GO  TO  9o2 
10  RETURN 

End 


SUBROUTINE  CAL  ( NS  .» ER.  SET#  RMSEH  ) 

COMMON  MS  >PAN»  B . ALPHA  A ( 33 ) »A(3  -O  ) t C < 3 5 ) » 

1 G(3  )»X(3'#5  )» ERROR (33#  3) # E( 5 # ) » XSET (3  ) 

2 » ALPHA  I(3>)  » YT  AB ( 20 ) » 

3 ALPZL ( 20 ) » NY  * CLAM 

b t CMUT AB ( 20 ) »CFOCIN(20) t DELFTB( 20 } » DELJTB (20 ) »KWAKE 

6 # EPST  AB ( 20 ) » LAM. LAMHL# BETA » SGAM  CLOT AB (20 ) 

7 »CN0uIil(20)  #DELNTB(20) 

DIMENSION  ER  ‘SET (33) 

REAL  LAM, LAMHL 
PI=3.14  S92oU 
RAD=57.29t> 

SBETA=SIM( BETA/RAD) 

MCP=NS  • 

M=2*NS  -1 

IF (NS  PAN. EG. 2)  NCP=M 
THETAS (LAM-BETA) /HAD 
FLIP-1.0 
DO  10  I=1»NCP 
Y=ABS('5*.S»C0S(PI*I/  (M+l ) 

IF ( I . GT . NS  ) THE  T AS ( LAM+BET  A ) /HAD 
IF' I.GT.NSS)  FLIP=-1.0 
IF  (I. EG.  MS  ) THETA=L-\M/NAO 
IE  ( I.EQ.NS  ) FLIPS').  • 

CTHETA=LOS( ThEtA) 

E PS=G I RC ( Y » r T AB » E PS  T Ab  * NY  * 1 ) / R AU 

CMUSGIPC ( Yf  YTALi.  LMUT  AB*  MY  1 ) *C0S»  EPS) / ( CThET  A* , 2 ) 
PLAM=C0S(  LAM-LAMHL)  / -<AO)  / COS  ( LAMHL/RAu ) 

CFO-.:  = GIPC  ( Y » YTAB.  CFOClNfUY  # 1 ) 

CNOC  = GIPC(Y» < T AH» CUOL IN. NY  < l ) 

UELFSSOiRC  (Yt  n AM,  DELE  TU.NY  . I) 

UELFzATAMl  TA;i(  UELPS/RAU)  *PLAv)  .RAO 
UELM=GIRC(  Y.  YTAB.UtLMT  >.NY»  l ) 

OELNsAT  •!!(  TAIPOELM/PAD)  •/Cl.AM)  *HAJ 
DEL  jSsf.  I ■>(.  ( Y • YTAM.DEL0Ti.NY.  I ) 

O'.LJSAT  N { T A M ( OEO S ♦ J‘  L JS ) / RAO ) *^LA^ ! m'AQ-  ELF 

CL  Os  GIRL  ( Y » YlAti.CLOTAM.rjY.  I ) 

i\E  z ( XSST!  I)  *M,\J»FLIM*v.r  TA4MGAM  )/CTMfc.TA) 

alp  I - ( alpha  ( i j ♦klIw*oLTa»Sg-\w  i/crm.  taw-h 
CAL  MACKE  T ( * ’AkE,  lFCC  .DLL'-  * CuPv  »Utl.M»  AE  » A(.M  I « l 1 « OSLO* 
1 OF  L A(  F » v O'.  -*) 

C 

C CL  IS  < IW-1/L-‘TI0M  LIFT  ONLY 

r 

7)  U.  = v L0»2.  i*f>I  * i l.i/CO  .ElAlf  ) 

c»*  I ) sCl  f t C The  r a * 2)  ‘C(  1i*o.*»/j 

1>  CONI  I Hi  ‘E 

CAL  WAT  VYJA.O  .LPM.A  I , NCP.  NLP.  1 . 3 » » '•  i 

00  20  I Si »MCP 

► P SET  (I  i = ( itSET  5 I ? } -ALPHA  t I?  \ ♦ Al/PMA  1 5 I ) *V/.29  M 
li  CONTINUE 
S’  . ■! 

DO  30  I isl.iiCP 


'3* 


I 


301  SUM=SUM+ER  'SET ( 1 1) * 2 
XNCP-NCP 

RMSER  ^(SQRKSUM/XNCP)  >*rad 
40  i RETURN 
END 
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SUBROUT  I f IE  HACKET  ( KWAKE  • CFOC  » DElF  » CNOC r DELNt  ALPE  » ALPI » CMU  * DEL  J» 
1 DELALF # COCP) 

EQUIVALENT  JET-FLAP  PROGRAM  USING  CLOSED  FORM  SOLUTION 
OF  JACOBS  EQUATION. 

INPUT* 

CF  IS  PHYSICAL  T.E.  FLAP  CHORD  (POSITIVE) 

DELF  IS  T.E. FLAP  ANGLE » V/.K.T  CHORD r JEGRrS  S NOSE  DOWN 
ALPE  IS  EFFECTIVE  INCIDENCE  RELATIVE  TO  CHORDLINE  C 
CMU  IS  FLAP  BLOW  IMG  MOMENTUM  COEF-  IC ! ENT t BASED  ON 
PHYSICAL  CHORD  C 
C DELJ  IS  ANGLE  DETW.O  N FLAP  CHORD  AND  JET 

DIMENSION  A ( 4 ) » PS II ( 4 ) # F ( 4 ) 

ALFDEL( X) =2.0/3. 14 ! S9*  ( SQRT( X* ( 1 .U-X)  + AS  I N( SQRT ( X ) 

PI=3.14B9 

KAD=PI/1S0. 

IF  (CMU)  2b‘J»CO 

2 CO  P=I.() 

CFPOCP=CFOC 
GO  TO  2(J 

3 IF  (CFOC)  2S0»4»rj 

4 CFOC=0.  1 

5 PI2=PI*PI 
PI32PI2*PI 
I T — U 

PS  i 1 ( 1 ) -U , 

XP-t'I/2 . 0 
PS  1 1 (3) =PI 
XDELF=RAD*DELF 
X0ELN=RAQ*  ;ELN 
XDELJ=RAD* Dfc LJ 
XALPE=RA0*ALP£ 

ETA=1 ,-OFOC 
ETAf-PCiiOC 
I X ALPI -RAD* ALP  I 
12  PS  I L i 2 1 - XP 

SlNXPSSlNt  XF>) 

C .’(RITE  (o»>9  >)  ITtXP.  (PS II  (F  ) »<*.  si.  3) 

’ FORMAT  ( 14  » ( 2 X * F l U . b ) 
tOSXP-CQM  XP) 

IT=  St  + 1 

T r E T A * ( i . +L USXP5-1 , 

P/ SAC OS ( T ) 

T 1 An* ( 1 . 4LQSXP)-1. 

POSACOStU) 

Hl):-I  xOCU  i~x  As,  PE  ) 

At  1 i-Pj 
F (2 ) SXDFLM 
A(2)=?»N 
F(3)=XDFLF 

A t J 5 

CHE- .«:»■'  ' 1 i *At  S)  * ! XR-SINXPi 
DO  no  ?P 2.3 

ASR-SlNIU.S* { A N 5 ♦ X R ) 
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ASM-SIN(0.bMA(N)-XP)  ) 

CHEC  K=C HEC K+  F { N ) * ( A ( N ) *XP+XP*S IN  < A ( N ) = -A ( N > *S  JNXP 
$+2.*ASM*ASP*AL0G( ASM/ASP) ) 

.1  0 CONTINUE 

CHECk=CHECK+XDELJ*XP*XP 

ANSI=CHECK/(1,0+COSXP) 

ANSJ=0 ,25*PI  t:CMU*  (SIN{  XDELF+XDELJ+XALPE) 

+ -SIN( (KWAKE-1. )*XALPI) ) 

RES=ABS ( ANSI ) -ABS( ANSJ) 

TOL=.0  =1 

IF  ( ABS  ( RES  ) -TOL ) 1 0 , L : ,15 
15  IF ( RES ) 40,10  ),30 
30  XP=(PSI1(  HPSI1  (2)  )/2.0 
PSI1(3)=PSI1(2) 

GO  TO  50 

40  XP=(PSJI (2)+PSIl (3) )/2.0 
PSIK1  )=PSI1  (2) 

50  CONTINUE 

IF ( ( IT-50) ) 12  t 12  1 98 
98  WRITE (6,10  ) RES, TOL* ALPE, XP 

1000  FORMAT ( * SOLUTION  FOR  EQUIV.  JET  FLAP  HAS  AN  ERR=*,F8.5, 
1 * T0I=*,F7.5,*  XSTAR=*,F10,5» * PSI1= 1 > F10 .5 ) 

10U  CONTINUE 

COCP- ( 1 . O + COSXP) /2  # 0 
CFPOCP-{1«0-T)/2«0 
CNPOCP- ( 1 . 0-T1 ) /2  = 0 
EMFOCP-  ( .1 . O-COSXP)  /2 . 0 
2 0 1 5 ADF=ALFOEL ( CFPOCP) 

ADEMF=ALFDEL( EMFOCP) 

ADN=ALFDEL(CNPOCP) 

DELALF-  ( 1 . O/COCP)  * ( ADN-1 . 0 ) *XDELN+ADF*Xl)ELF+ ADEMFtXDEL  J) 
250  RETURN 
END 
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SUBROUTINE  FIELD  ( Nil'  NPTS' B#  ALPHAB*  BETASS*  SWEEP#  GAMA# 

1 PHI  I # GAM AO V # XPB # YPB  # ZPB > MP2 ) 

PROGRAM  FOR  ESTIMATION  OF  THE  VORTEX- INDUCED  SUBSONIC-FLOW  FIELD 
BENEATH  SWEPT  AND  UNSWEPT  WINGS. 

REFERENCE  NACA  REPORT  1327  BY  WILLIAM  J.  ALFORD#  JR. 

PROGRAM  BY  A.  E.  HOLMES  OF  DEPT  72-74*  -- *FEB.  1967*m- 

UPDATED  ON  JAN  IS  1971 

UPDATED  FOR  SUBROUTINE  FOR  JET  FLAP  hIGH  LIFT  PROGRAM 
23  OCTOBER  1973.  AEH 

CHANGED  TO  HANDLE  ASYMMETRICAL  LOADINGS  25  OCT  73.  AEH 
VELOCITIES  INSIDE  CIRCLE  OF  RADIUS  S ARE  ZEROED  OUT.  29  OCT  73 
LAST  UPDATE  27  NOV  73. 

DIMENSION  XP ( 20 ) » YP(20)#  ZP(20)»  X(50)#  Y(50)# 

1 Z (50 ) » GAM0V2 ( 50 ) 

2 » PHI  I (35 ) # GAMAOV ( 3 ) » XP3 ( 20 ) # YPB (20 ) » ZPB ( 20 ) 

PI=3. 14159 
RAD-180 . /PI 

TGAM=SIN( GAMA/RAD) /COS (GAMA/RAD) 

tlam=sin(Sweep/rad)/cos(swe:.p/rad) 

NVORTS=50 
VORSPA=B/50. 

S-VORSPA/2 . 0 
BETArl.O 

DO  40  J=1»NV0RTS 
Y ( J) =S-( B/2 . 0 ) + ( VORSPA* ( J-l ) ) 

ABSY-ABS ( Y ( J) ) 

Y W PH I -AC OS  (2,0*Y  (J)/B) 

GAM0V2 { J ) =G IRC ( YWPHI » PHI  I # GAMAOV . MP2 » 2 ) 

X(  J)-A!3SY*TLAM 
Z(J)=0.  i 
69o9  FORMAT ( ) 

40  CONTINUE 

SALF=SIN( ALPHA B/RAD) 

CALF-COS ( ALPHAB/RAO ) 

SBETA=SIN(  BETAS’ ./RAD) 

CHET A- CoS( BETAS'  /RAD) 

DO  50  I — 1 » NPTS 

YP( t)=XPB(I)*SBETA+YPB( I)*CBETA 

XP< I)=ZPB( I ) *SALF+ !XPB( I ) *CBETA-YPU( I)*SUETA) *calf 
7P( I ) -ZPB ( I)*CAl^-XPB( I)*SALF 
50  CONTINUE 

WRITE  <0*11 1 ) 

1 FORMAT  (/  # 5 X ► * X POINT  Y POINT  7 POINT  *# 

1 » EPSILON  SIGMA*#/) 

DO  20  1=1# NPTS 

SUMGFW=U 
SUMGFV=0 
SUMGFU-0 

DO  10  J=1 ♦ NVORTS 

DELX- ( XP ( I ) -X ( J) ) /BETA 
QELYR=YP( I )-Y ( J) 

DELZ- ZPt I)-Z ( J) 

S-VORSPA/2. U 
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DN^vJ=5QRT(  ( DELX/S)  *-2+  (DELZ/S)  **2  ) 

DNSQMR=SQRT  ( (DNSQ**2)+  ( ( DELYR/S ) -1 ) * 2) 

SSJELR=(  ' DELYR/S )~1  )/DNSQMR 
C -DELR=DNSQ/DNSQMR 

DNSOPR=SQRT ( ( DNSQ*  2 ) + ( ' DELYR/S ) + 1 ) * *2 ) 

SN8ETR=( ( DELYR/S) +1 i/DNSQPR 

CSBETR=DNSQ/DNSQPR 

SNGAMAr ( DELX/S) /DNSQ 

CSGAMA= ( DELZ/S) /DNSQ 

IF ( DNSQ. LT • 1 » 0 ) GO  TO  62 

FWRPL= ( SNGAMA/DNSQ ) * ( SNBETR-SNDELR ) 

GO  TO  64 
62  FwRPL-0  . 

64  IF ( ! DNSQPR.LT. 1.0). OR. (DNSQMR.LT.l-0)  GO  TO  78 
FWRPL-F'VRPL- 

1 ( ' DELYR/S) -1 ) / ( DELZ/S  5 * 2+(  DELYR/S)-1)*  2 ) > * ( 1+SNGAMA* 

3 CSDELR  ) + ( •'  t DELYR/S  ) 4-1  )/(  ! DELZ/S  ) * ; 2 + 

4 t ( DELYR/S) +1) **2)  ) * (1+SNGAMA*CSBETR ) 

FVKPL=( ( DELZ/S)/ ( DELZ/S)*  2+(  DELYR/S ) +1 )*  2 ) Ml  + SNGAMA 

2 *CSBETR)-< i DELZ/S) /(  DELZ/S)*  2+ 

3 ( ( DELYR/S ) “1 ) * ♦ 2 ) ) * ( 1+SNGAMA*CS0ELR ) 

GO  TO  80 

78  FVRPL-U , 

80  IF ( DNSQ. LT. 1.0)  GO  TO  8 

FUWPL=( CSGAMA/OMSQ) * ( SUBETk-SNUELR) 

GO  TO  90 
8 FURPL=0.  > 

90  SUMGF,-/=SUM6F-.V+(GAM0V2(  J)*FwRPL) 

SUMGFV=SUMGFV+ ( GAM0V2 ( J) *FVRPL ) 

SUMGFU=SUMGFL'+(GAM0V2  ( J!  *FURPL) 

1.0  i CONTINUE 

.vAOV-SUMGF  ,V/  ( 4 ♦ R I fS ) 

VA0/=SUMGFV/(4*PI* j) 

UA0w=SUM6Fii/{4*Pl  »S) 

FPSL0N3  ( ATAf-K  a AQV/  ( 1 fUAQV  ! *87.3 

SIGMA=-{ AT  \M( VAQV/( 1+UAOV)  : *67,3 
OLOO-  ( 1 +UAOV ) * *2  * t V AO  / i * 2*aA0V* 

.VMITE  (6.1  0)  X PU  ( I ) * YPU  ( I ) * ZPB(l).  EPSLOU#  SIGMA 
110  FORMAT  (hX.RFIO.S) 

20  continue 

RF  TURN 
E ND 
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